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ABSTRACT 


A research  program  was  conducted  to  develop  a validated  numerical 
model  for  the  prediction  of  radiation  emitted  in  the  infrared  wavelength 
spectral  band  from  turbojet  and  rocket  exhaust  plumes.  Two  existing 
gas  dynamics  models  (LAPP  and  REP3)  were  employed  to  calculate  the  flow 
field  properties  (temperature  and  species  concentrations)  required  as 
input  for  the  radiation  computer  codes.  A new  finite  difference  computer 
code,  which  was  developed  as  part  of  this  project,  and  includes  the 
calculation  of  imbedded  shock  waves  in  the  exhaust  plume  along  with 
finite-rate  chemistry  and  turbulent  mixing,  was  only  partially  successful, 
A refined  version  of.  the  MIRADCOM  Radiation  Model  was  developed  and  was 
used  successfully  to  predict  CC^  and  ^0  radiation  emitted  from  several 
cases  for  which  experimental  data  were  available. 


1.  INTRODUCTION 


During  the  past  few  years,  considerable  effort  has  been  expended  to 
design  infrared  seeker  systems  for  ground-to-air  and  air-to-air  missile 
guidance  systems.  Since  the  next  generation  of  seeker  systems  must  possess 
head-on  capability,  infrared  radiation  from  the  missile  exhaust  plume  is 
an  important  component  of  the  overall  infrared  signature  of  a missile  system. 
In  order  to  characterize  completely  the  missile  seeker-target  interaction, 
it  is  necessary  to  develop  infrared  radiation  predictive  models  which  have 
the  capability  of  predicting  all  the  pertinent  physical  phenomena  associated 
with  the  plume  infrared  signature. 

Infrared  radiation  emitted  from  the  hot  gases  in  an  exhaust  plume  is 
characterized  by  the  types  of  radiating  species  and  their  number  densities 
and  the  environmental  conditions  (temperature  and  pressure) . All  radiative 
predictive  models  require  a knowledge  of  these  parameters,  so  that  the  first 
step  in  developing  a model  is  to  develop  or  to  obtain  a gas  dynamics  model 
which  yields  the  spatial  distribution  of  radiating  specie  concentrations, 
temperature,  and  pressure. 

Unfortunately,  the  exhaust  plumes  of  turbojet  aircraft  and  rockets  are 
characterized  by  the  existence  of  a complex  shock  wave  structure  and  the 
presence  of  a turbulent  viscous  mixing  layer  between  the  inner  hot  jet  flow 
and  the  external  freestream  flow.  The  inner  flow  usually  consists  of  a 
high-temperature,  high-velocity,  chemically-reacting  gas  produced  by  the 
combustion  of  a hydrocarbon  fuel  while  the  external  flow  usually  is  low- 
temperature,  low-velocity  air  that  flows  around  the  exterior  of  the  aircraft 
or  missile  body. 


Even  though  considerable  progress  had  been  made  in  developing  workable 
gas  dynamics  computational  models ,[1-17]  at  the  time  this  project  was  initiated, 
there  was  no  available  model  which  adequately  predicted  the  mixed  supersonic 
and  subsonic  flow  regions  present  in  the  central  core  of  the  exhaust  plumes. 
This  report  describes  the  results  of  a research  program  which  was  initiated 
to  develop  a finite  difference  computer  code  which  would  include  the  calcula- 
tion of  imbedded  shock  waves  and  Mach  discs  along  with  finite-rate  chemistry 
and  turbulent  mixing  phenomena. 

During  the  course  of  this  project,  a gas  dynamics  computer  code  developed 
by  Dr.  Brian  Spalding  for  the  Rocket  Propulsion  Establishment  in  the  United 
Kingdom  was  made  available  to  the  U.S.  Army  MIRADCOM  as  part  of  the  TTCP 
program  with  Great  Britain.  This  program  was  modified  by  the  Principal 
Investigator  to  operate  on  the  MIRADCOM  computer.  The  results  of  a few 
limited  calculations  are  described  in  this  report. 

The  MIRADCOM  Radiation  Model  (developed  by  H.  Tracy  Jackson,  Jr.  and 
extended  by  K.  E.  Harwell)  was  refined  as  part  of  the  research  effort  described 
in  this  report.  Results  obtained  using  this  model  will  be  described  and  a 
comparison  made  between  the  theoretical  predictions  and  experimental  data. 
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2.  THEORETICAL  FOUNDATIONS 


2.1  Gas  Dynamics  Theoretical  Models 

A complete  review  of  the  theory  for  turbulent  exhaust  plumes  is  beyond 
the  scope  of  this  report.  A description  of  the  flow  theories  is  given  by 
Jackson,  Poslajko,  and  Harwell  in  Ref.  18, 

The  appropriate  equations  of  motion  for  the  gas  are  the  conservation 
equations: 

3 o — k 

Conservation  of  Mass:  g-~-  + V • pV  = 0 

Conservation  of  Momentum: 

P [fir  + (V  • V)  v]  = - Vp  + V • x 


- V . (pV)  + V • JT  • V - V . q 

After  supplementing  these  conservation  equations  with  an  Equation  of 

State 


Conservation  of  Energy: 


:>■  DE  „ 
p — = P 
P Dt 


»E 


+ v 


Dt 


e = e (T,p) 


and  constitutive  relations  relations  relating  the  shearing  stress  T to  strain 
of  the  fluid  element  and  ultimately  to  the  velocity  gradients  in  the  flowj 
a complete  system  of  flow  equations  is  obtained.  For  example,  the  laminar  shear 
stress  can  be  expressed  as 


iJ 


'a! i*»!i 


X V • V q.  + 2p 
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The  theoretical  approach  to  the  development  of  the  turbulent  shear  stress 
model  is  not  described  here,  but  the  model  employed  was  the  eddy  diffusivity 
model  proposed  by  Donaldson  and  Gray  [ 19]  for  the  compressible  free  mixing 
of  a primary  jet  with  quiescent  air: 


T = — p u'v' 


„-LU 
K p s 


3u 

3r 


where  K is  the  mixing  rate  factor,  L is  a typical  local  scale  length.  Pis  the 

local  density,  and  U is  a typical  velocity  difference.  These  parameters  are 

s 

defined  in  Refs.  [ 12  and  19] . 

It  had  been  intended  to  incorporate  the  Turbulent  Kinetic  Energy  (TKE) 
turbulence  model  approach  into  the  computer  model.  However,  due  to  problems 
with  the  basic  computer  code,  this  was  not  accomplished.  In  this  model  [see 
Ref. 11  ],  an  additional  conservation  equation  for  the  turbulent  kinetic  energy 
is  added  to  the  previous  set  of  conservation  equations.  This  equation  for  an 
axisymmetric  flow  becomes 


3k  , - 3k 

P;U3X  + Pv37 


+ pe. 


3k 
't  3r 


3u 

3r 


V2 


- a 


2 Jl, 


where 


k = 


(u-)2  + (vT)2 


, 3 u 

h - Vp  3F 


al p k 


and 
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et 
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a^  and  are  Universal  empirical  relationships  nnd  is  a length  scale  for 
the  TKE.  The  reader  is  referred  to  Harsha  [ 20]  for  a detailed  formulation 
of  this  model.  Suffice  it  to  state  that  the  addition  of  the  TKE  equation  in 
the  matrix  formulation  of  the  computer  model  is  straightforward  and  can  be 
accomplished  without  major  difficulty. 

The  chemical  reaction  relations  employed  in  the  gas  dynamics  model 
was  identical  to  those  employed  in  the  Low  Altitude  Plume  Program  (LAPP) 
developed  for  the  U.S.  Air  Force  Rocket  Propulsion  Laboratory  [ 5 ].  The 
reader  is  referred  to  the  LAPP  report  [ 5 ] for  a detailed  description  of 
the  formalism. 

2.1.1  Gas  Dynamics  Finite  Difference  Computer  Model 

A finite  difference  numerical  scheme  was  used  to  approximate  the  gas 
dynamics  conservation  equations  described  in  the  previous  section. 

A method  developed  by  Daniel  Matuska  of  the  USAF  Weapons  Laboratory 
for  computing  hydrodynamic  properties  of  fluids  was  modified  to  include  the 
viscous  stress  tensor. 

The  finite  differencing  technique  used  was  the  standard  one  in  which  the 
derivative  of  the  function  F is  defined  as  [F(X  + H)  - F(X)]/H  where  H is 
finite. 

The  degree  of  success  in  using  the  finite  difference  technique  to  solve 
the  differential  equations  depends  on  the  "stability,"  "convergence"  and 
"compatibility"  of  the  method.  The  method  is  stable  if  any  bounded  starting 
procedure  yields  a uniformly  bounded  solution  of  the  difference  equation  as 
H goes  to  zero.  A method  is  convergent  if  the  values  of  Y(N),  the  solution 
of  the  finite  difference  expression  approximating  the  difference  equation, 
tend  to  the  values  of  Y(T)  of  the  exact  solution.  The  method  is  compatible 
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if  the  truncation  error  goes  to  zero  as  H goes  to  zero.  The  difference  between 
the  true  solution  of  the  differential  equation  and  the  finite  difference 
equation  method  is  influenced  by  the  step  size  H and  may  require  very  small 
values  of  H to  be  stable  and  convergent. 

One  of  the  major  advantages  of  the  finite  difference  approach  is  that 
the  problem  is  reduced  to  an  initial  value  problem  with  well  defined  boundary 
conditions  which  results  in  a unique  time-dependent  solution.  The  major 
advantage  in  terms  of  the  present  problem  of  interest,  is  that  the  mixed 
subsonic/transonic/supersonic  flow  problem  can  be  solved  without  using  a 
patching  solution. 

The  main  disadvantages  of  the  approach  are  that  the  complete  boundary  must 
be  defined  as  a function  of  time,  that  small  steps  may  be  required  which 
results  in  a fine  computational  mesh  which  in  turn  requires  a large  computer 
memory  if  realistic  resolution  is  desired,  and  that  computational  times  may 
become  prohibitive. 

The  calculation  procedure  used  in  the  present  program  is  to  treat  the 
problem  in  two  phases.  In  the  first  phase,  the  "Lagrangian  terms"  (coordinate 
system  fixed  to  the  fluid  element)  are  integrated.  In  the  second  phase,  the 
convective  terms  are  integrated  (mass  momentum  and  energy  is  transported  across 
boundaries) . The  details  of  the  calculation  procedure  are  described  in 
Ref.  [21]. 


-6- 


2.2  Infrared  Radiation  Theoretical  Model 


A numerical  infrared  radiation  band  model  was  developed  to  predict 
the  infrared  spectral  and  spatial  radiation  intensity  distributions  in  exhaust 
plumes.  The  model  essentially  extends  the  models  of  Jackson  [22,23],  which 
were  developed  for  CC>2  radiation  from  aircraft  plumes,  to  include  water  vapor 
radiation  and  to  include  the  capability  of  treating  higher  temperatures  than 
the  Jackson  code. 

The  statistical  band  model  and  spectral  parameters  used  in  the  code 

were  developed  for  application  to  rocket  plume  studies  [ 24  ] . The  transmissivity 

in  the  model  becomes 
K 

tr  = exp  [-P/(l  + P2/4Q fy 

where  P = / kdu 

Q = / ~ du 

k(w,T)  = absorption  coefficient 
u = optical  depth 
Y = average  line  width,  and 
d(«,T)  = average  line  spacing 

The  absorption  coefficients  k and  reciprocal  line  spacings,  1/d,  ate  tabulated 
as  a function  of  frequency,  0)  and  temperature,  T in  Ref.  [24], 

The  average  line  width  was  computed  using  the  relation 

r " [120.1  PH2o/T]  + t.09<P„2o  + PS2)  + .04  P02  + .12  PC02!  [273/T]15 

2.2.1  Infrared  Radiation  Computer  Model 

The  infrared  radiation  computational  model  combined  a water  vapor 
emissivity  formalism  developed  by  Reardon  [25]  with  the  plume  integration 
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scheme  developed  by  Jackson  [ 23  ] for  carbon  diokide  emission  from  aircraft 
plumes. 

The  computation  of  radiance  from  a jet  exhaust  plume  was  accomplished 
by  dividing  the  plume  into  sufficiently  small  incremental  blocks  so  that 
average  conditions  could  be  assigned  and  the  blocks  treated  as  homogeneous. 

The  approach  used  in  the  geometric  division  of  the  plume  is  described  in 
detail  by  Jackson  [23  ].  The  plume  is  subdivided  into  slabs  parallel  to  the 
line  of  sight.  Each  slab  is  then  cut  into  horizontal  strips  which  lie  along 
the  line-of-sight  vector. 

The  radiation  from  each  of  these  strips  is  computed  by  integration  using 
the  modified  Curtis-Godson  approximation  [24]. 

Integration  is  accomplished  over  the  length  of  the  strip  to  obtain  the 
parameters  P and  Q as  shown  earlier, 

P = / kdu 

n r kY  , 

Q = / — j-  du 

„ p- 

The  transmittance  T = exp  [~P/(1  + p2/4Q)2]  is  then  computed. 

R 

Radiance  is  then  calculated  by  integrating  over  the  strip  emissivity. 

For  a homogeneous  block,  Radiance  = (Blackbody)  (emissivity)  so  that  the 

inhomogeneous  strip  radiance  is 

STRIP  RADIANCE  = / (Blackbody)  d£ 

Since  de  = - dx  the  computational  algorithm 
K 

STRIP  RADIANCE  = - £ (Blackbody) .Atr  . 

i * 

The  algorithms  used  for  calculating  the  radiation  from  inhomogeneous  strips 
were  tested  by  applying  them  to  homogeneous  strips  and  comparing  with  direct 
radiation  calculations.  The  differences  in  calculated  values  were  less  than 
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one  percent. 

A listing  of  the  infrared  computer  program  is  given  in  Appendix  I. . 

A list  of  the  input  requirements  is  given  in  Appendix  II.  Appendix  III  includes 
a listing  of  a sample  case  with  output. 


JL DISCUSSION  OF  RESULTS 


This  section  contains  a discussion  of  the  results  obtained  as  part  of 
the  contract  effort.  The  procedure  used  to  validate  the  infrared  radiation 
computer  code  is  described  in  Section  3.1.  Section  3.2  presents  the  results 
of  the  infrared  radiation  code  predictions  for  a turbojet  flight  test  for 
which  there  are  experimental  data.  The  prediction  of  the  infrared  radiation 
emitted  from  a small  kerosene/gaseous  oxygen  rocket  engine  exhaust  plume  is 
presented  in  Section  3.3.  A discussion  of  the  progress  in  developing  the  Gas 
Dynamics  Finite  Difference  Computer  Model  is  described  in  Section  3.4.  Some 
preliminary  results  obtained  using  the  REP3  Computer  Code  are  described  in 
Section  3.5. 

3.1  Validation  of  the  UTSI  Infrared  Radiation  Code 

In  order  to  validate  the  Infrared  Radiation  Code,  the  spectral  radiance 
of  water  vapor  and  carbon  dioxide  mixtures  experimental  data  are  available 
[Ref s. 24, 27-29] . Sukanek  and  Davis  [26]  performed  a similar 
assessment  of  the  NASA  band  model  formulation  developed  by  Reardon  [ 25  ] . 

It  was  decided  that  several  selected  cases  used  by  Sukanek  and  Davis  would  be 
used  for  the  evaluation  since  this  would  provide  a comparison  of  the  present 
model  calculations  with  the  experimental  data  and  with  the  predictions  of  the 
NASA/Reardon  Radiation  Code.  The  cases  selected  for  the  comparison  are  given 
in  Table  3.1. 

The  spectral  radiance  of  60-cm  thick  slabs  of  hot  H^O  vapor  and  CO2  were 
computed  for  the  2.7-p  and  4.3-y  wavelength  bands.  Figures  3. 1-3. 8 present  the 
results  of  the  computations. 

The  comparison  between  the  measured  and  predicted  values  of  the  spectral 
radiance  in  the  2. 7- U band  for  an  isothermal  slab  of  water  vapor  is  presented 
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TABLE  3.1  CASES  CONSIDERED  FOR  COMPARISON 


Isothermal 

Hot 

H20 

Isothermal 

Hot 

CO? 

Isothermal 

Hot 

CO  2 

Isothermal 

Mixture 

Isothermal 

Mixture 

Nonisothermal 

H?0 

Nonisothermal 

C02 

Particle  Pressure  (mmHg) 


Temperature  Profile 


Nonisothermal  Mixture 


760 


28  675 


382/537/723/953/1128/1160/990/751/558/389 


386/528/719/953/1130/1160/975/737/541/387 


378/537/723/958/1127/1158/990/752/555/383 


NOTE:  The  temperature  of  the  gas  was  either  homogenous,  in  which  case  one  temperature  is  reported,  or  was 
inhomogeneous,  in  which  case  the  temperatures  are  given  at  the  center  of  each  6 cm  wide  cell. 
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in  Fig.  3.1.  At  the  center  of  the  spectral  band 'the  computed  values  are 
slightly  lower  than  the  measured  ones.  The  NASA  code  predicts  slightly 
higher  values  than  the  UTSI  code,  but  in  general  the  agreement  is  fairly 
good  between  the  calculated  and  measured  values.  The  calculated  values  are 
within  20%  of  the  measured  ones. 

The  measured  and  predicted  values  of  spectral  radiance  in  the  2,7-11 
and  4,3-y  wavelength  bands  for  an  isothermal  slab  of  CO2  are  presented  in 


Figs.  3,2  and  3.3.  In  contract  to  the  U^O  predictions,  the  predictions  of 


the  spectral  randiance  for  CC^  do  not  agree  very  well  with  the  measurements. 
Errors  of  30  to  50  percent  are.  observed.  Similar  disagreement  was  found  by 
Lindquist  and  others  [27-29]  who  attributed  to  the  probable  cause  of  the  dis-^ 
agreement  to  the  lack  of  knowledge  of  the  CO2  broadening  parameters  at  low 
pressures.  Ludwig  and  others  [ 24]  show  better  agreement  between  the  pre- 


dictions and  measurements  at  higher  pressures  (P^  3 1-2  atmospheres) , 


Figures  3.4  and  3.5  present  the  results  for  the  spectral  radiance  of  an 


isothermal  mixture  of  CO2  and  1^0  in  the  4,3  and  2,7-y  bands,  respectively, 


The  calculated  values  are  higher  than  the  measurements  in  the  center  of  the 
spectral  bands.  Both  codes  yield  results  which  are  in  fair  agreement  with 
the  measured  values. 

Figure  3,6  presents  the  results  for  a nonisothermal  slab  of  water  vapor, 
There  is  fairly  good  agreement  between  the  predicted  and  measured  values,  but 
the  present  code  appears  to  predict  values  which  are  below  the  measured  and 
calculated  values  using  the  NASA  code. 

Figure  3.7  presents  the  comparison  between  the  predicted  and  measured 
values  of  spectral  radiance  in  the  2.7-y  spectral  band  of  a nonisothermal 
slab  of  002*  The  agreement  between  the  predicted  and  measured  values  is 
quite  good. 
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Fig.  3.4  SPECTRAL  RADIANCE  IN  THE  4.3-ym  BAND  FOR  AN  ISOTHERMAL  MIXTURE  OF 


CO2  and  ^0. 
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Fig.  3.5  SPECTRAL  RADIANCE  IN  THE  2.7-jlm  BAND  FOR  AN  ISOTHERMAL  MIXTURE  OF 
CO.  and  H.O. 


SPECTRAL  IRRADIANCE 


Figure  3.8  presents  the  results  for  a nonis'othermal  mixture  of  C02  and 

H^O.  There  is  fairly  good  agreement  over  the  entire  2.7-y  wavelength  band. 

In  summary,  the  UTSI  radiation  code  appears  to  yield  results  which 

agree  fairly  well  with  the  NASA  code  even  though  there  are  some  differences. 

With  the  exception  of  the  low  pressure  CC>2  cases,  the  predicted  values  are 
within  ± 20%  of  the  measured  values. 


WAVE  NUMBER  (cm-1) 


Fig.  3.8  SPECTRAL  RADIANCE  IN  THE  2.7-ym  BAND  FOR  A 
NONISOTHERMAL  MIXTURE  OF  C02  AND  H20. 


3.2  Infrared  Radiation  Predictions  for  Turbojet  Flight  Test 


The  UTSI  Infrared  Radiation  Code  was  utilized  to  calculate  the  infrared 
radiation  emitted  in  the  wavelength  interval  from  4 to  5 microns  from  a 
turbojet  aircraft  exhaust  for  which  experimental  data  are  available. 

H.  Tracy  Jackson,  Jr.  at  MIRADCOM  provided  the  data  given  in  Table  3.2. 

The  Low  Altitude  Plume  Program  (LAPP)  was  exercised  for  the  input  data 
given  in  Table  3.2.  The  first  portion  of  the  LAPP  code  output  for  this  case, 
labeled  "FLTl-PASS  7,"  is  included  as  Figs.  3.9-3.11.  The  output  data  from 
the  LAPP  code  was  then  used  as  input  data  for  the  Infrared  Radiation  Code. 

Refer  to  Fig.  3.12  for  a typical  set  of  input  data. 

The  results  of  the  computations  are  presented  in  Figs.  3.13-3.15.  The 
CC>2  contribution  to  the  radiation  per  length  of  plume  (Watts/SR/cm)  is  given 
in  Fig.  3.13  while  the  H^O  contribution  is  given  in  Fig.  3.14.  The  I^O  radiation 
is  very  small  compared  to  the  CC^  radiation. 

The  spatial  distribution  of  the  radiation  is  given  in  Fig.  3.15.  Also, 
shorn  on  the  figure  is  the  centerline  temperature  distribution.  The  total 
radiation  emitted  from  the  plume  in  the  4.0  to  5.0  micron  wavelength  interval 
was  165.7  Watts/SR.  This  radiation  is  at  the  plume  and  is  uncorrected  for 
detector  response  characteristics  or  for  atmospheric  attenuation  between  the 
plume  and  the  detector. 

It  had  been  anticipated  that  the  Jackson  MIRADCOM  Radiation  Code  would  be 
used  to  calculate  the  radiation  for  the  same  test  flight  so  that  a comparison 
could  be  made  between  the  two  codes.  Unfortunately,  the  results  from  the 
MIRADCOM  Code  did  not  become  available  prior  to  the  completion  of  the  project. 
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TABLE  3.2 

# 
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Initial  Conditions  for  LAPP  Computer  Program 
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Fig.  3,.  10  Sample  Output  of  LAPP  Computer  Program 
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Fig.  3.12  Typical  Set  of  Input  Data  for  Radiation  Code 
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Fig.  *.3.13  ' CO2  Contribution  to  Infrared  Radiation  per  Length  of  Plume 
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3 . 3 Infrared  Radiation  Predictions  for  a Small 
Kerosene/Gaseous  Oxygen  Rocket  Engine 

Two  gas  dynamics  and  the  UTSI  infrared  radiation  computational  models 
were  used  to  predict  the  physical  properties  in  hot  rocket  exhaust  plumes. 
Experimental  data  were  obtained  in  the  exhaust  plume  of  a small  kerosene/gaseous 
oxygen  rocket  engine  using  a three-dimensional  laser  Doppler  velocimeter  system 
to  measure  the  velocity  distribution.  The  spatial  distribution  of  radiated 
energy  in  the  4-5  y wavelength  band  was  measured  using  an  infrared  radiometer. 
Detailed  comparisons  between  the  computational  and  experimental  results  are 
described  in  this  section. 

з. 3.1  Description  of  Experimental  Equipment 

Small  Laboratory  Rocket  Engine:  The  exhaust  plumes  studied  in  these 
measurements  were  produced  by  a small  kerosene/gaseous  oxygen  rocket  engine  [ 30  ] 
operating  in  the  NASA  Laser  Doppler  Velocimeter  Facility  [ 31  ] . Engine  design 
and  performance  parameters  are  listed  in  Table  3.3.  Data  were  taken  at  a 
variety  of  oxidizer/fuel  (0/F)  ratios  and  at  two  chamber  pressures.  The  ex- 
perimental data  used  to  compare  with  the  theoretical  calculations  were  obtained 
at  a chamber  pressure  of  11.1  atm.  and  an  0/F  ratio  of  2.25. 

NASA  Laser  Doppler  Velocimeter  Facility:  The  NASA  Marshall  three-dimensional 
laser  Doppler  velocimeter  [31]  was  used  to  measure  the  velocity  distribution 
in  the  exhaust  plume  and  mixing  layer  produced  by  the  inner  hot  rocket  exhaust 
jet  mixing  with  an  outer  subsonic  airflow.  A schematic  diagram  of  the  LDV 
experimental  arrangement  showing  the  installation  of  the  small  rocket  engine 
is  shown  in  Fig.  3.16.  As  described  in  Ref.  [18]  the  simultaneous  measurement 
of  the  three  Doppler  velocity  components  yielded  the  three  velocity  components 

и, v,  and  w. 
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TABLE  3.3 


ENGINE  DESIGN  AND  PERFORMANCE  PARAMETERS 


Fuel 

Kerosene 

Cxidizer 

Gaseous  oxygen 

Purge 

Nitrogen  100  psig 

Coolant 

Water  at  3 gal/min 

Chamber  pressure 

50  and  148  psig 

Fuel  flow  rate 

50  pslg 

1.38  g/sec 

148  psig 

3.80  g/sec 

Oxygen  flow  rate 

50  psig 

3.05  g/sec 

. 148  psig 

8.70  g/sec 

Throat  diameter 

0.200  in. 

Exit  diameter 

0,310  in, 

Throat  area 

0.0314  in.2 

Exit  area 

0.0755  in.2 

Area  ratio  (A^/A^) 

2.40 

0/F  ratio 

2.25 

Chamber  diameter 

0.780  in. 

Chamber  length 

1.620  in. 

* 

L 

25.20  in. 

C at  148  psig 

6045  fps 
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In  order  to  simulate  external  airflow  aroutid  the  small  rocket  entine, 
an  outer  flow  nozzle  was  used  to  produce  an  external  flow  velocity  of  61  m/sec. 

A schematic  diagram  of  the  rocket  engine  and  outer  flow  nozzle  is  presented 
in  Fig.  3.17. 

Infrared  Radiation  Diagnostic  Equipment : The  primary  instrument  used  to 
obtain  the  radiation  intensity  data  was  an  Electro-Optical  Industries  Model  470 
Radiometer.  The  data  were  obtained  using  a fixed  filter,  which  transmitted 
radiation  in  the  4-5  y band. 

A large  metal  shield  with  a 1-mm  aperture  was  moved  in  the  axial  direction 
to  obtain  the  variation  of  radiant  intensity  along  the  exhaust  plume  centerline. 
A second  set  of  measurements  was  made  using  a 1-mm  vertical  slot  in  the  large 
metal  shield.  The  slot  length  was  greater  than  the  diameter  of  the  exhaust 
plume.  The  experimental  setup  is  shown  in  Fig.  3.18. 

3.3.2  Typical  Experimental  Results 

Velocity  Distribution  Data:  Experimental  data  for  mean  flow  velocities 
(u,v,w),  turbulent  intensities  [(u')2,  (v')2,  (w’)2],  and  turbulent  velocity 
correlations  were  obtained  in  the  radial  direction  at  various  axial  locations 
downstream  from  the  nozzle  exit  plane.  Data  were  obtained  at  X/D  locations 
of  2.4,  4.8,  8.4,  11.3,  and  14.2,  where  X is  the  distance  measured  from  the 
nozzle  exit  plane  and  D is  the  diameter  of  the  exit  nozzle.  Typical  data 
obtained  at  X/D  = 11.3  for  the  mean  axial  flow  velocity  are  presented  in 
Fig.  3.19.  Other  data  are  described  in  Ref.  [18]. 

Infrared  Radiation  Data:  The  infrared  radiation  data  have  been  described 
previously  [30  ].  Typical  data  obtained  at  an  0/F  ratio  of  2.25  and  an  ex- 

ternal flow  velocity  of  210  fps  are  presented  in  Figs.  3.20  and  3.21  for  the 
1-mm  aperture  and  slot,  respectively.  The  centerline  radiation  intensity 
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Fig.  3.18  Exhaust  plume  test  arrangement  in  NASA  LDV  facility. 


-35- 


U MEAN  VELOCITY  (ft/sec) 


ODATA  FROM  POSITIVE  Z POSITIONS 
©DATA  FROM  NEGATIVE  Z POSITIONS 


W,Z  N INTERNAL  JET 
EXHAUST 

RADIUS  (Re!  = 0.155  in. 

NOTE:  CHAMBER  PRESSURE 
AT  148  psig 


© © © © © 


0 1 2 3 4 5 6 7 

Z/R„ 
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Fig.  3,20  Spatial  distribution  of  radiation  along 
axis  of  plume  for  aperture  diameter  = 1 
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Fig.  3.21  Spatial  distribution  of  radiation  along 
axis  of  plume  for  1-mm  slot  width;  0/F  = 
velocity  = 210  fps. 


distribution  in  Fig.  3.20  clearly  shows  the  presence  and  influence  of  shocks 
and/or  Mach  disks  in  the  exhaust  plume.  The  effect  of  afterburning  in  the 
mixing  layers  on  the  infrared  signature  is  shown  clearly  in  Fig.  3.21  where 
the  shock- induced  radiation  peaks  are  superimposed  on  an  axially  increasing 
radiation  intensity.  The  data  presented  in  Figs.  3.20  and  3.21  were  obtained 
in  the  near  field  of  the  exhaust  plume  (X/D  less  than  12)  and  indicate  that 
the  peak  radiation  intensity  has  not  been  reached.  Typically,  the  axial  dis- 
tribution of  radiation  is  similar  to  that  shown  in  Fig.  3.22.  Even  though 
the  data  presented  in  Fig.  3.22  are  for  a lower  mass  flow  and  chamber  pressure, 
the  axial  variation  is  similar  to  that  obtained  for  the  high-chamber-pressure 
cases.  The  data  presented  here  are  only  samples  selected  to  indicate  the 
measured  trends  in  preparation  for  a discussion  of  the  comparison  between  theory 
and  experiment. 

3.3.3  Comparison  of  Calculated  and  Measured  Results 

The  LAPP  and  TKE  computer  codes  were  exercised  for  the  rocket  performance 
parameters  listed  in  Table  3.3.  The  NASA  Lewis  Chemical  Equilibrium  Computer 
Program  was  used  to  calculate  the  exit  plane  gas  properties.  Because  of  the 
short  expansion  nozzle  length,  the  chemical  species  were  assumed  to  be  frozen 
at  the  throat  values.  These  exit  plane  conditions  then  were  used  as  initial 
conditions  for  the  LAPP  and  TKE  computer  codes. 

The  results  predicted  by  the  LAPP  and  TKE  computer  codes  are  presented 
in  Figs.  3.23-3.28.  The  variation  of  the  centerline  velocity  with  distance 
from  the  nozzle  exit  plane  is  presented  in  Fig.  3.23.  The  LAPP  code  predicts 
a constant-velocity  core  having  a length  approximately  equal  to  X/D  = 5,  which 
is  followed  by  an  exponential  velocity  decay  due  to  turbulent  mixing  predicted 
by  the  Donaldson-Gray  viscosity  model.  The  TKE  (REP3)  code  predicts  a sharp 
'decrease  in  velocity  at  approximately  an  S/D  value  of  one.  There  are  also 
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Hg.  3,23  Calculated  axial  velocity  distribution 
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several  other  oscillations  in  velocity,  as  shown  by  the  jagged  nature  of  the 
curve.  A particularly  large  decrease  in  velocity  occurs  neat  an  X/D  value  of 
8.  Also  noted  in  Fig.  3.23  are  the  locations  of  the  measured  radiation  peaks 
shown  previously  in  Fig.  3.20.  The  X/D  locations  for  the  LDV  measurements 
also  are  indicated  on  Fig.  3.23. 

The  predicted  variations  of  temperature  at  the  jet  centerline  with  axial 
-distance  are  given  in  Fig.  3.24.  The  LAPP  code  predicts  a small  temperature 
rise  due  to  afterburning  in  the  exhaust  plume.  The  TKE  code  predicts  a sharp 
rise  in  temperature  at  an  X/D  value  of  one,  which  then  is  followed  by  a series 
of  oscillations  in  temperature.  The  predicted  oscillations  do  not  agree  with 
the  measured  locations  of  radiation  peaks  (denoted  by  arrows  on  the  Figs,). 

It  is  quite  obvious  that  neither  of  the  computer  codes  does  an  adequate  job  of 
predicting  the  spatial  variation  of  velocity  and  temperature  produced  by  shock 
waves  and  Mach  disks  in  the  near  field  of  the  exhaust  plume.  This  is  not 
surprising  since  LAPP  has  no  provision  for  treating  shocks. 

Figures  3.25-3.38  present  the  variation  of  the  axial  component  of  the 
velocity  with  radial  distance  from  the  jet  centerline.  As  shown  in  Fig,  3,25 
for  an  X/D  value  of  2.4,  the  predicted  results  (curves)  do  not  agree  with  the 
experimental  measurements  (symbols).  Because  of  the  limited  number  of  particles 
in  the  hot  core  of  the  jet,  the  NASA  LDV  system  was  not  able  to  make  measurements 
closer  to  the  jet  centerline  than  about  one  exit  nozzle  radius.  The  variation 
of  the  radial  velocity  component  also  is  shown.  The  TKE  code  predicts  an  in- 
creasing value  of  v with  radius,  whereas  the  measured  value  decreases  with  radial 
distance.  The  external  flow  is  an  axial  flow  at  61  m/sec. 

The  lack  of  agreement  between  theory  and  experiment  also  is  indicated  in 
Fig,  3,26  for  the  velocity  profile  at  an  X/D  value  of  8.0.  The  LAPP  and  TKE 
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Fig.  3.25  Velocity  profiles  at  X/D  = 2.4 
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Fig.  3.26  Velocity  profiles  at  X/D  = 8.0 


calculated  results  are  similar  but  are  in  excess  of  the  experimental  values. 
Similar  results  are  presented  in  Figs.  3.27  and  3.28  for  X/D  values  of  11.3  and 
14.2. 

The  LDV  data  indicate  that  the  exhaust  plume  is  much  narrower  than  pre- 
dicted by  the  LAPP  and  TKE  gas  dynamics  models . Since  both  the  LAPP  and  TKE 
models  assumed  that  turbulent  mixing  began  at  the  nozzle  exit  plane  and  ignored 
the  large  recirculation  (base  flow)  region,  it  is  expected  that  turbulent 
mixing  occurs  too  fast  compared  to  the  physical  situation  where  turbulent 
mixing  begins  at  the  end  of  the  recirculation  zone,  which  is  several  nozzle 
diameters  downstream  from  the  exit.  It  is  known  that  the  base  flow  region 
can  have  a pronounced  effect  on  exhaust  plume  spread  and  initial  turbulent 
mixing  layer  growth.  Because  of  time  and  budget  constraints,  we  were  unable 
to  use  a base  flow  model  to  assess  its  effect  on  the  calculated  temperature 
and  velocity  profiles. 

The  infrared  radiation  band  model  program  was  used  to  calculate  the 
variation  of  infrared  radiation  along  the  exhaust  plume  for  the  temperature 
and  species  concentration  profiles  predicted  by  the  LAPP  and  TKE  gas  dynamics 
computer  codes.  The  computed  results  are  shown  in  Fig.  3.29.  The  lack  of 
agreement  between  theory  and  experiments  was  not  unexpected,  since  the  calculated 
temperature  distributions  did  not  yield  the  spatial  structure  due  to  shock 
waves  that  are  evident  in  the  measured  infrared  data  and  are  visible  in  photo- 
graphs of  the  exhaust  plume. 

It  is  not  clear  what  causes  the  difference  between  the  REP3  and  LAPP 
calculations.  The  LAPP  code  plus  the  radiation  code  yield  results  that  are 
within  a factor  of  1.5  of  the  measured  results  while  the  factor  is  approximately 
2.5  for  the  REP3  plus  radiation  code  results. 
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Fig.  3.29  Comparison  of  calculated  and 
measured  IR  radiation. 
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3.4  Gas  Dynamics  Finite  Difference  Computer  Model 


A Gas  Dynamics  Finite  Difference  Computer  Model  was  developed  as  part 
of  the  contract  effort.  Numerous  difficulties  and  problems  were  encountered 
during  the  development  process,  but  the  program  was  operational  at  the  conclusion 
of  the  effort.  Unfortunately,  the  present  version  of  the  code  is  not  entirely 
satisfactory  due  to  the  fact  that  the  computational  times  are  large.  Approxi- 
mately one  hour  on  the  IMB  360  computer  system  was  required  to  march  the 
solution  360  cycles  (which  corresponds  to  a flow  time  of  0.2  msec.).  While 
this  time  would  be  approximately  6-10  minutes  on  the  CDC  660  computer,  it  is 
still  felt  to  be  quite  excessive  since  the  test  flow  being  analyzed  was  a relatively 
simple  expanding  flow.  In  its  present  form,  the  finite  difference  code 
successfully  calculates  shock  waves  on  wedges  and  blunt  bodies  and  is  moderately 
successful  in  calculating  the  shock-Mach  disc  structure  in  a cold  gas  nozzled 
expansion. 

Several  typical  nozzle  expansion  test  cases  were  computed  to  assess 
whether  the  numerical  results  were  reasonable.  These  results  are  described 
in  the  next  section. 

3.4.1  Overexpanded  Nozzle  Exhaust 

In  order  to  assess  the  FD  model’s  ability  to  handle  strong  shock  waves 
which  result  in  the  exhaust  plume  from  an  overexpanded  nozzle,  the  flow  having 
the  initial  characteristics  given  in  Fig.  3.30. 

The  computed  Mach  number  variation  along  the  jet  centerline  is  given  in 
Fig.  3.31.  The  Mach  number  decreased  to  a subsonic  value  in  a very  short 
distance  downstream  from  the  exit  plane  and  remained  subsonic  for  80-cm  in  the 
flow  direction. 

An  analytical  estimate  of  the  flow  conditions  was  made  by  assuming  that 
an  oblique  shock  wave  would  be  formed  at  the  nozzle  exit  and  would  raise  the 
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pressure  behind  the  shock  wave  to  the  ambient  pressure  of  1.0  atm.  The  oblique 
shock  angle  was  then  determined  to  be  approximately  76.6°  as  shown  in  Fig.  3.32. 
The  shock  wave  would  intersect  the  centerline  axis  at  a distance  of  0.2323-cm 
downstream  of  the  nozzle  exit  plane.  Gas  dynamics  calculations  yielded  a 
value  of  0.94  for  the  Mach  number  in  region  (2)  behind  the  shock  wave.  This 
calculated  value  is  indicated  on  Fig.  3.31  and  is  in  good  agreement  with  the 
finite  difference  code  result. 

As  shown  in  Fig.  3.32,  the  shock  wave  would  reflect  from  the  centerline 
and  would  process  the  flow  back  parallel  to  the  centerline  axis.  In  doing  so, 
the  pressure  in  region  (3)  would  increase  while  the  Mach  number  M^  would  be 
less  than  M2.  This  agrees  qualitatively  with  the  FD  code  results. 

At  the  intersection  of  the  reflected  shock  with  the  free  boundary,  an 
expansion  wave  would  be  generated  to  reduce  the  pressure  at  the  edge  of  the 
plume  back  to  one  atm.  The  expansion  fan  would  accelerate  and  expand  the  flow 
in  regions  (4-6)  until  a new  shock  pattern  would  be  formed  between  regions 
(7  and  8).  Then  the  entire  pattern  would  repeat. 

While  the  FD  code  calculations  appear  to  be  quantitatively  correct,  there 
were  no  available  numerical  data  to  compare  with  the  FD  code  results  in  the 
subsonic  flow  region.  In  order  to  determine  whether  the  FD  code  predictions 
are  reasonable  in  an  expanding  flow,  another  test  case  was  designed  in  which 
a supersonic  flow  would  expand  from  a nozzle.  The  results  for  this  flow  are 
described  in  the  next  section. 

3.4.2  Underexpanded  Nozzle  Exhaust 

As  stated  above,  a test  case  involving  the  supersonic  expansion  of  air 
from  a nozzle  was  employed  so  that  the  FD  code  predictions  could  be  compared 
with  the  results  obtained  using  the  Method  of  Characteristics. 


Fig.  3.32  ASSUMED  SHOCK  PATTERN  FOR  HAND  CALCULATION. 
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Fig.  3.33  COMPARISON  OF  METHOD  OF  CHARACTERISTICS  AND  FINITE  DIFFERENCE  MODEL  RESULTS 


The  method  of  Characteristic  technique  used  was  Temple's  Method  as  taught 
by  the  Principal  Investigator  in  a first-year  graduate  course  in  gas  dynamics. 

The  characteristics  net  produced  in  the  X-Y  physical  plane  is  shown  in  the 
lower  portion  of  Fig.  3.33.  A comparison  of  the  MOC  results  (using  5 characteristics 
to  start  the  calculation)  and  the  Finite  Difference  Code  results  is  presented 
in  the  upper  portion  of  Fig.  3.33.  The  results  are  in  good  agreement.  The  lack 
of  agreement  at  large  distances  from  the  nozzle  exit  plane  is  to  be  expected 
since  the  hand-calculated  MOC  solution  used  only  five  characteristics  in  the 
initial  fan  and  is  therefore  expected  to  be  more  approximate  than  the  Finite 
Difference  Computer  Code  which  uses  a more  complete  set  of  differential  equations. 

Even  though  we  were  encouraged  that  the  computer  code  appears  to  yield 
reasonable  numerical  results  for  an  expanding  plume,  we  are  discouraged  by  the 
fact  that  the  computer  program  required  approximately  one  hour  on  the  IBM  360 
computer  system  to  march  the  solution  0.3  msec,  (or  360  cycles).  While  this 
time  would  be  approximately  6-10  minutes  on  the  CEC  6600,  it  still  is  felt  to 
be  quite  excessive  in  view  of  the  simple  flow  calculated. 

3.4.3  Other  Test  Cases 

In  earlier  work  [21],  the  finite  difference  code  was  used  to  calculate 
the  viscous  flow  between  two  parallel  walls  with  one  of  the  walls  moving 
relative  to  the  other.  This  is  the  so-called  classical  Couette  Flow  problem 
whose  theoretical  solution  is  well  known  [32],  The  computed  numerical  results 
agreed  well  with  the  theoretical  model  results.  In  addition,  the  computer 
solution  converged  in  10  to  20  iteration  cycles  and  the  total  computational 
time  on  the  IBM  370  computer  was  approximately  30  seconds. 

The  code  also  successfully  predicted  the  flow  over  a two-dimensional 
wedge  in  which  a strong  bow  shock  wave  is  formed. 
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In  both  of  these  cases,  we  have  solid  boundaries  which  restrict  the 
flow  and  appear  to  assist  in  the  convergence  of  the  solution  to  a stable  value. 
3.4.4  Chemical  Reactions  and  Turbulent  Mixing  Computations 

Due  to  problems  with  the  basic  finite  difference  computer  code,  the 
chemical  reaction  and  turbulent  mixing  portions  of  the  code  were  not  fully 
operational  at  the  end  of  the  program. 

3.5  Results  Obtained  Using  the  KEF  3 TKE  Computer  Model 

The  REP3  TKE  Computer  Model  described  earlier  was  revised  for  use  on 
the  CDC  computer  at  MIRADCOM.  The  test  case  supplied  by  Dr.  David  Jensen  of 
the  Rocket  Propulsion  Establishment  of  Great  Britain  was  successfully  run  oh 
the  MIRADCOM  computer. 

After  the  REP3^  program  was  operational,  it  was  utilized  to  calculate 
the  flow  field  from  a small  kerosene/gaseous  oxygen  engine  similar  to  the  one 
used  as  the  test  case  by  Jensen.  The  results  of  these  calculations  are  pre- 
sented in  Figs.  3.23-3.28.  As  discussed  in  Section  3.3,  the  REP  3 predictions 
do  not  agree  with  the  experimental  data.  Attempts  to  locate  any  problems  with 
the  calculations  or  with  the  experimental  data  were  unsuccessful. 

One  computer  run  was  devoted  to  rerunning  the  case  described  in  the 
previous  paragraph.  The  program  was  rerun  using  a much  smaller  grid  spacing. 
The  results  are  shown  in  Fig.  3.34.  The  top  curve  shows  the  oscillations  in 
temperature  while  the  lower  curve  indicates  the  results  for  the  velocity. 

Also  shown  on  the  graph  are  the  locations  of  the  radiation  minima  (two  experi- 
mental runs  are  shown:  Runs  65  and  85).  As  can  be  seen,  the  first  three 
locations  of  the  calculated  temperature  minima  agree  with  the  experimental 
locations.  Even  though  the  fourth  through  seventh  minima  do  not  agree,  the 


results  are  reasonable. 
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However,  the  behavior  of  the  calculated  profiles  for  X greater  than  90-mm 
are  disturbing  in  that  there  are  a series  of  temperature  and  velocity  oscillation 
which  do  not  agree  with  any  experimental  or  physical  evidence.  Furthermore, 
the  predicted  sharp  drop  (rise)  in  temperature  (velocity)  at  145-rom  and  the 
sharp  rise  (drop)  in  temperature  (velocity)  at  165-mm  do  not  correspond  to 
observed  physical  phenomena. 

Clearly,  further  efforts  are  needed  to  ascertain  whether  or  not  the  REP3 
Program  is  correct.  Unfortunately,  there  was  no  available  time  or  funds  to 
consider  this  further. 

It  was  also  discouraging  that  over  30  minutes  of  computer  time  on  the 
CDC  6600  was  required  to  obtain  the  small  step  size  required  to  resolve  the 


near  field  shock  waves. 


4.  SUMMARY  AND  CONCLUSIONS 


A research  project  was  conducted  which  resulted  in  the  development  of  an 
extended  version  of  the  MIRADCOM  (Jackson)  Infrared  Radiation  Computer  Model, 

The  extended  computer  code  was  validated  by  comparing  the  predicted 
spectral  radiance  of  ER^O  and  CO2  in  the  2.7-ym  and  4.3-ym  wavelength  bands 
with  experimental  data  and  with  predictions  of  the  NASA  (Reardon)  Radiation 
Computer  Code.  The  predicted  results  agree  fairly  well  with  the  NASA  code 
even  though  there  were  some  slight  differences.  With  the  exception  of  a few 
results  for  low  pressure  CC^,  the  predicted  values  were  within  ± 20%  of  the 
measured  values. 

The  infrared  computer  code  was  used  to  predict  the  infrared  radiation 
emitted  in  the  wavelength  interval  4 to  5 microns  from  the  exhaust  of  a full 
scale  turbojet  aircraft.  The  predicted  results  appeared  to  be  reasonable 
when  compared  with  the  available  test  data. 

The  infrared  radiation  band  model  was  used  to  calculate  the  variation 
of  infrared  radiation  along  the  exhaust  plume  of  a small  kerosene/gaseous 
oxygen  rocket  engine.  The  temperature  and  species  concentration  profiles 
were  predicted  using  the  LAPP  and  REP3  computer  codes.  The  agreement  between 
the  theoretical  calculations  and  the  experimental  measurements  was  poor.  The 
lack  of  agreement  was  due  to  the  fact  that  the  calculated  temperature  distributions 
did  not  exhibit  the  spatial  structure  due  to  shock  waves  that  are  evident  in 
the  measured  infrared  data  and  are  visible  in  photographs  of  the  exhaust  plume. 
There  appears  to  be  short  comings  in  the  ability  of  both  the  LAPP  and  REP3 
codes  to  predict  the  correct  spatial  distributions  of  temperature  and  specie 
concentrations. 
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The  Gas  Dynamics  Finite  Difference  Computer  Model  developed  as  part  of 
this  project  was  only  partially  completed  and  was  able  to  produce  reasonable 
numerical  results  at  the  expense  of  long  computer  runs  and  high  computer  memory 
requirements.  The  code  successfully  calculated  the  flow  field  for  flow 
geometries  having  solid  boundaries,  but  had  problems  with  the  free  boundaries 
encountered  in  nozzle  exhaust  plumes.  The  results  obtained  for  the  near 
field  of  the  nozzle  appeared  to  be  correct,  but  the  downstream  results  were 
not  validated. 

The  preliminary  calculations  obtained  using  the  REP3  and  TKE  computer 
codes  did  not  agree  with  the  experimental  data. 

Due  to  the  partial  success  of  the  finite  difference  and  REP3  computer 
codes,  further  development  and  assessment  are  required  before  these  codes 
are  fully  operational  in  a validated  model. 
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APPENDIX  I 


LISTING  OF  IR  COMPUTER  PROGRAM 


IIOIXII  III  III  IIXIQTI 

//HARWE  JOB  ( 369)1 ,51 1 1 8 . REM3TE3 . R- 1 92 . T- 35 . h- 1 0 ) , NGUYEN  : 

/ / STEP  l EXEC  EORTGCLG,FnRTREG=192K ,GOREG=I92K 
//FORT  .SYSPRINT  DD  DUMMY 

//fORT.SYSIN  DD  * 

COMMON / PLOT  1 /WAVEI , W A VEF , CENTR , TOTRAD 
COMM  0 1'J  /RASP/RASP 

COMMON  /BAND?/  XKK  ( 2 5 4 , 7 ) , DD  I N V ( 2 5 4 , 7 ) , UK  ( 2 5 4 ) , A KU(  254)  J._ _____  

DIMENSION  1 6 (3), T 8 ( 1 0 ) , 17(6),  12(10), 13(6), 14(4), 15(10), II C 5) 

COMMON  / R P M 1 / N W A V , N A V N 0 ( 2 5 4 ) , W A V [,  T ( 2 5 4 ) 

COMMON/CHECK/TREAD  _ _ 

COMMON  /INDATA/  X ( 4 0 ) , R C 4 0 , 3 0 ) , T (40 , 30  ) , PC02  ( 40,30) , 

$ PH  20(40, 30)  , P N 2 ( 4 0 , 30)  ,P02(40, 30)  , PCO(4(),  30)  , PH  2 ( 40, 30) 

C 0 M M.  0 N _ / C T R 0 1,  / _T  I Tf . E ( 1 2 )_,  K 0 U N Tj.N  X S ,_N  0 R U N 

1 COMMON  /OTH/  PP  , R.j , NRP  ( 40  ) 

COMMON  /SWITCH/  KSWIT  ,IXXJG 

c _ _ _ _ ■ - _ 


DIMENSION  R ADSTL.  ( 254  ) , RAD  XT  (25  4 ) , TRLAST ( 254  ) 000130 

DIMENSION  CR(IOO)  , C X ( 1 0 0 ) , CL (100)  000140 

DIMENSION  R4DXt,(254)  ,RRN(  100)  __ 

DIMENSION  STR AD (40 ) 

DIMENSION  RAD (2 54) ,RADX(254) ,TRNS(254) ,P(254) ,0(254) ,RADSTR(254)  000230 

DIMENSION  WAVE 5(254)  _000?40 

DIMENSION  RCA  DC (10  0)  ' ’ " 000260 


DIMENSION  SR  AD (2 54) ,SRADSR(254) 

INTEGER  I H 20(2)  ,ICD2(2)  , I WRT(?) 

REAL,  HRADXT(40)  , C R A D X T ( 4 0 j " ' 

REAL  HCTOTX(40)  ,HCRPX1’(40) 

__  __  REAL  HR  ADX  ( 254 ) ,CRADX(254)  ,SRADX(254)  __  

EQUIVALENCE  ( HR  ADX  ( 1)  , R ADX  ( I ) ) , (CR  AOX  ( 1 ) , RAD  ( T)  j',  

S ( SRADX ( 1 ) , TRNS ( 1 ) ) 

EQUIVALENCE  ( HCTnTX ( J ) , CR ( 1 ) ) , ( HCRDXT ( 1 ) ,CX(1)) 

c ""  ’ ■“ . " " ' 

DATA  1 6 N / 1 0 / 

DATA  I 8 M / 3 7 / 

” DATA  I7N/23/  ' ' 

DATA  1 2 N / 3 3 / 

D A T A I 3 N / 2 5 / 

DATA  I4M/1P/  : r ' • '■ 

DATA  I 5 N / 4 9 / 

DATA  I1N/20/ 

DATA  ASTER , «CALC/t H* /i  00*1  H 7 ” '000300 

DATA  I DISC/ 8/ 

DATA  I STOP / 0 / 

'data  iDvSCT/o/ ’ " ' " 

. DATA  J PASS/O/ 

DATA  IELG/0/ 

ISNIT  = 2 ' " " 

kswit  = 2 
J S W I T = 2 

LSWIT  = 1 ' 

L R E A i)  = 0 
KPZQ=0 
NZU=1 
KN  IQ  = 0 
IXXJGsO 

C ' 

C SENSE  SWITCH  1 SET- - YES , PR  1 NT  FLOW  FIELD  DATA 

C SENSE  SWITCH  ? SC-T--YES,  CQ2  PADI ATION  — NO,  H20  PADIATIOM 

C SENSE  SWITCH  3 SET  — YES,  AERDCHKM  INPUT  TAPK  — HU,  motor  DATA  input  tape 

C SENSE  SvITCM  4 SET-- YES, H/O  AND  CO?  RADIATION  COMBINED 

C 
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n n o o i ! o i ! | no  onnn 


ror«AD=o.o 


000590 


SET  PRINT  LABEL 


DETERM INE_  MODE„OF  TMPIIT_  TAPE,.--.  MOTOR  . DATA  OR  AEROCHEM 

172  CONTINUE 
ENl)  FILE  9 

REWIND  9 * 7' 

END  FIDE  8 

REWIND  8 _ 

IF(KPZQ.K0. 1)  GO  TO  104  ' ‘ 

GO  TO  (101 , 102)  , JS-VIT 

102  CONTINUE  _ _■  _ _ 

call'  reader ~ ' ~ ~ "" 

KPZQ=1 
GO  TO  104 

1 01  CONTINUE-  “ ' " -----  - 

104  CONTINUE 

NXSiKOUNT  

DO  111  1=1, NXS 

IF(NRP(I)  .GT.25)  GO  TO  1 1 1 

NRP  1 =NRP( I ) + l 
I SET=NRP ( I ) 

’ DO  112  K K = N R P 1 ,25  __ " 

■'■R(I,KK)=«(i',ISET) ~ 

T( I , KK ) =T ( I , ISET) 

PC02  (I  ,KK  ) =PC02  ( I , ISET)  _ ___ _ 

PR20  (I  , KK ) =PH20( I , I SET ) " 

PN2 ( X , KK ) =PN  2(1, ISET) 

P02(I  ,KX)=P02(T  ,ISET)  - __ 

112  C 0 U TIM U E 
111  CONTINUE 

“READ "WAVELENGTHS  AMD  COR RES POM DING  BAND  PARAMETERS  ~ 

READ  IM  H20  BAND  DATA  FIRST  FOLLOWED  BY  C02  BAND  DATA 

47  CONTINUE 

IFCLSWIl’  .EO.t  .AMD.  ISrOP.EO,  3.  AND.IFLG.LT.  2)  GO  TO  21 
IF ( I STOP , EG , 3 ) GO  TO  2 3 
21  CONTINUE 

WRITE  (6,9911)  000330 

READ  (5,9912)  MWAV  l , NMAV2  _ 000340 

'WRITE  (6,1  0035)  MWAV1  , MWAV2 
10035  FORMAT ( 1 OH  NWAV1  ,?(1X,I5)) 

9912  F0RMATC2I  10)  . 

N W A V = 1 + ( ( N W A V 2 - N 4 A V 1 ) / 2 5 ) 000 360 

DO  .1  1 = 1 , N W A V 000380 

IF(KNTO.EO.O)  GO  TO  9957  • 

GO  TO  9967 

9957  READ  (5,9922  ) WA VMO( 1 ) , ( XKK ( I , J ) , J = 1 , 7 ) 

9922  FORMAT  (F9.l,7E0.2)  000400 

G 7 I’O  9977 

9967  READ  (6,9922  ) W A V;i  0 (I)  , ( X KK  ( 1 , J ) , J = 1 , 7 ) 

9977  CONTINUE 

WRITE  (6  , 100  37  ) »AV  IO( f ) , (XKK ( J , J) , J=|  , 71 
10037  FORMAT  ( 1 Uii  WAV  ,F9.1,7E9.2) 

(F(K-IT').EO.O)  GO  TO  R 9 8 7 
GO  TO  '< q 9 7 

9°R7  READ  (S,‘iOOt)  v A V " 0 ( T ) , ( D D I N V ( J , J ) , U = 1 , 7 ) 

9001  FORM  AT ( F6, 1 ,E9.2,6E1 0.2) 
ll  i'll  99  17 

9997  READ  (6,9010)  •*  A Vn  0 ( T ) , f o'  D I D V ( I , J ) j J = l , 7 ) 

9010  FORM/vf  (!■•/;,.  1 , K9.?.  oEIO.  2) 
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9947  CONTINUE 

WRITE  (6,10037)  WAV»lfJ(I),(DDINVCr#J),J  = l,7)  * 

WAVLT(  U = 10000. /WAVMD( I)  000420 

9910  FORMAT  ( 1 X , RE 1 5 . 7 ) 000450 

1 WAVL5(I)=WAVLT(I)**5  000460 

4AVeisWAV.40(l  ) _ ' 000470 

W A V E F=  W A V N O ( rl  VI A V ) , ’ ”■  000480 

WRITE  (6,9911)  000490 

9911  FORMAT  (1H1)  _ __  ....  000500 

00  9913  1 = 1 , N W A V 000510 

00  9913  0 = 1 ,7  000520 

_IF  ( ODIN V ( 1,0)  .OF.  0,0)  ODIHV( I , J)  = l .0E-20  J000530 

IF( XKK(  I ,0)  , IjE.O,  ) XKK(I,J)=1  .0E-2Q 

9913  CONTINUE  000540 

9000. FORM AT^ F9 . 4 , 7E9 , 3 ) 


PRINT  BAND  PARAMETERS  IF  OPTIONED  SELECTED 


IF  (LSW.IT  . EQ.t , AMO.  ISTDP . EG  . 0 . AND . IFl.AG . I.T  . 1 ) CALL  SAVE1 
I F ( I S W I T .EO.l)  CALL  BPPRT 

23  CONTINUE 

I F ( L S W I T ,EQ.l.ANO.ISrOP.EQ.3.AND.IFLG.GE.2)  CALL  SWAP 
IFLG  = IFLG+ 1 


DETERMINE  INPUT  DATA  PRINT  OPTIONS 

SUBROUTINE  PRMD  WRITES  OUT  INPUT  FLOW  FIELD  DATA 
IFdSWIT  .EO.l)  CALL  PRMD 


ALPHA= ASPECT  ANGLE  IN  DEGREES (0  DEG=NOSE-ON ) 

IST0P=1 ,WJLL  TERMINATE  PROGRAM  AFTER  SINGLE  ALPHA  CALCULATION 
IST0P=2, RETURNS  TO  THIS  POINT  FOR  ANOTHER  ALPHA 
I STOPs 3 , RETURNS  TO  STMT  172  TO  READ  NEW  SET  FLOW  FIELD  D AT  1 
SENSE  SWITCH  4 WILL  OVER  RIDE  ANY  ISTOP  OPTION 

103  READ (5 , 1 0 1 0 ) ALPHA , I S TOP 
1010  FOR M AT (FI  0.5,15) 

ALP=ALPUA/57 . 295780 
CC=COS( ALP) 

SC=SIN( ALP) 

PI  = 3.  141  5926 

TRADK  = 0.0  •' 

TO  TR  AD=0 , 0 
CENTER=Q . 0 

DO  185  I 1 = 1 , N S'  A V 

185  RAI)XT(  ( 0=0.0  

WRITE  (6,1401) 

DO  3000  1=1, NXS 
NXPEX=NRP(T) 

WRITE  (6,1401) 

140  1 FORMAT  (1HO 

WRITE  (*>,JJ04)  ( T T T L K ( J R R ) , J R t<  = 1 , 1 2 ) 

1 404  F )R  LA  rtLU,  1 244  , A2,  /) 

VRI  TE  (6,  1.400)  I ,X(  O , ALPHA 

1 400  FORMAT  ( IX,  I GUSTATION  4 0 . = , 1 4 , 5 X , I 7 H A X I A l>  1)1. St.  ( C P ) = , F 1 0 . 3 , 5 X , 
*1  44  ASPECT  A ’•!  G I , <■'  (DEG)  = ,F5. 1 ,/) 

r r aox  = o . o 

1X9=1 

Rt)  = f<(  r (11,25) 

IF  (Rn  , L t1  . 0.0)  »n  = 1 . OH-20 

X )=<(  .IX") + .0  1 

I F(  ALPHA  .EO  .180.0)  X'1  = X(  MXS) 

L p » = i/.PiU 
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001520 

001530 

001540 

001550 

001560 

001570 

001580 

001690 

001600 

001610 

001620 

001630 

001640 

001650 


001670 

001680 

001690 

001700 

001710 

001720 

001730 

001740 

001750 

001760 

001780 
001800 
001910 
n 0 l 8 2 0 
001430 
00 1 940 
001 850 


♦♦♦♦change  delz=rii/zpez  ro  delz=i.  for  the  si^ap*** 

OELZ=i.O 

DELZ  = Ril/ZPEX 

*******  SET  MZU=NZPEX  FDR  TOTAL  PLUME  ***** 

ZU=R ASP/DHLZ 
NZU=MZPEX 


00  190  tt=l,NWAV 
SRADSR (II ) =0 .0 
190  RAOXC  II)=0.0 
STOT=0. 


001870 

0018R0 


C THIS  CARO  HAS  BEEN  CHANGED  AS 

C 00  2000  M=  1 , NZPEX 

DO  2000  M = 1,NZIJ  _ 

Z M = M 001910 

ZO=DKLZ*(ZM-l . ) 001920 

IE  ( ALPHA  . EQ  ,0,0)  R0=2.*R(1 ,25)+(ZM-l  . ) * ( R ( NXS  , 25) -2  , *R  (.Lr  2_5JJ./Z_PEXOO  l 9 30 
IF ( ALPHA. EQ . 180.0)  RO  = ( ZM- 1 . ) * ( R ( NXS , 25 ) ) / ZPEX  001940 

C 001950 

^♦♦♦♦♦♦♦♦S************^******************************* ***♦♦¥*♦♦♦♦♦ *****001960 
C CALCULATE  COORDINATES  X,R  ALONG  L INE-OE-SIGHT  AS  FUNCTION  OF  ITS  001970 

C LENGTH  001980 

C __ 001990 

* DELL=io7  


IFCALPHA.EO.O.O.OR. ALPHA. EG. 180.0)  DELL=X ( MXS) /98 . 
200  XL=0  . _.  .....  ...  


C 

C 

C 


C 

c 

c 

c 

c 

C 

C 


00  203  J=l,100 

N = iJ  . ..  . 

APR=XL*SC-SGRT (RO*RO-Z9*z6) 

CR(N)  -SORT  (ZO*ZO+  APR*  APR) 

CX(M)=XQ+XL*CC 
. C L ( N ) = X L 

lE(CXCM) . LT , 0 . ) GO  TO  204 
1F(CX(N> .GT.X(NXS) ) GO  TO  204 
IE (ALPHA. HQ, 0.0. OP, ALPHA. EG, 180.0)  GO  TO  203 
00  201  K = 1 , N X S 

K t = K ......  , 

IECCX(N)  .JjE.X(K)  ) GO  TO  202 

201  com  HUE 

202  IE(CR('1)  ,GT.R(KI,2S)  .AND.  N , GT . 10)  GO  TO  204 

203  XL=XL+DELL 

DELL  TOO  SMALL, INCREASE  AS  .FOLLOWS 


0 E L L = 0 E L L * 1 . 5 

GO  TO  200  _ _ ' 

204  CONTINUE 

1 E ( N , G E . 5 0 ) GO  TO  299 

DELL  TOO  LARGE, DECREASE  AS  FOLLOWS 

205  OELL=CL ( N ) /SO . 

PLU'lK  THICKNESS  -igGLICl  OLE  AT  THIS  POINT,  «OVF  Oh  TO  NEXT  POINT 
AFTER  SETTING  ALL  SPECTRAL  R ADT  6 IT  ON  VALUES  EQUAL  10  ZERO 

l S’ ( 0F|,I, , GT . .005)  G O TO  200 

210  c iurnoE 

00  21  1 1 1 = 1 ,<!!>’  a v 

RH'STLC  l I.  ) - 0 , <i 

211  RADS Tw( I I)  =0.0 
GO  TO  1990 


002000 
002010 
0.0202  0 

002030 
002040 
002050 
002080 
002070 
002080 
002090 
002100 
002110 
002120 
002130 
002140 
002150 
002170 
002  1 80 
002190 
002200 
002210 


002260 

002270 

002280 

002300 

002320 
002330 
0 0 2 3 (1 0 


002390 

002400 

002480 

002460 
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C*****+**********************+******************* ************** ********* 0 02470 
C CALCULATE  PRESSURE  AND  TEMPERATURE  ALONG  LINEfOE-SIGHT  AT  5ACH  X,R002490 

C 002490 

299  CONTINUE 

WRITE! 6,307) DELZ , DELL 

307  EORMAT(2(tX,n:tO„4)) . ..  __  

c 

DO  300  11=1, NWAV 

P(II)=0,0  _ 002510 

0 C 1 1 ) =0 . 0 

TRNS(II)sl.0  002530 

300  RADSl’R  ( II  ) =0 , 0 ___  _ 002540 

C 

DO  1000  0=1, M 002560 

J CHECK= J __  __  _____  _ 002570 

C 

DO  301  L = 2 , M X S 

If  (JCHECK.EO.99)  LCHECK=3  

HXl=L-l 
MX2  = L 

IF(CXCJ)  ,LE.X(L)  ) GO  TO  302 ' 

301  CONTINUE 

C 

302  CONTINUE  _ __  _ __ _____ 

DO  30  3 L = 2 ,25 

If (JCHECK.EO. 99)  LCHECK=3 

N R 2 = L _ _ __  __ 

If (CR(J) .LE.R(NX2,L) ) GO  TO  304 

303  CONTINUE 

"3 OT  CONTINUE  “ —— — — — - : 

DO  305  L=2 ,25 

NR1=L  ' . 

I F ( C R ( J ) . I,  E . R ( M X 1 , L ) ) GO  TO  306 

305  CONTINUE 

306  CONTINUE 
D R 1 = C R ( J ) - R ( U X 1 , M R 1 - 1 ) 

DR2  = R ( NX  1 , NR  1 ) "R  ( NX  1 , NR  1 - 1 ) _ , 

DR3  = CR ( J ) -R ( NX  2 , NR 2 - t ) 

D R 4 = R ( N X 2 , N R 2 ) - R ( r ! X 2 , N R 2 - 1 ) 

0X1=CX(,1)-X(NX1) 

DX2=X(NX2)-X(nxi) 


C 


C 


0 R 1 0 R 2 , DR  3 DO  4 , D X 1 D X 2 WERE  CREATED  TO  ELIMINATE  INDEFINITE  FORMS  002760 

Dill  DR  2 = 0.0  002770 

D R 3 0 R 4 = 0 . 0 _ . _ 002780 

0X10X2  = 0,0  ’ ~ ’ ' ~ 0027  90 

1 f ( 0 R 2 . N E . 0.01  DW1DR2  = 0R1/DR2  002H00 

If  (DR4.NE.  0,0)  QP3DR4=DR3XDR4  __  002919 

I f ( 0 X 2 .ME.  0.0)  0X10X2  = 0X1/0X2  ’ 0029  20 

T X 1 =T ( N X 1 , * I P 1 — J ) + ( V ( N X 1 , N R 1 ) -T ( U X 1 , NK 1 - 1 ) ) *D«1 DR 2 

T < ? = T ( N X 2 , NP2-1  ) + f T ( N X 2 , M R 2.  ) - T ( N X 2 , I R 2 - 1 ) ) * f R 3 D R 4 

C T = T X 1 + ( I X 2 - T X 1 ) * 0 X 1 D X 2 002.950 

RH2K  l=Pn2.n(  IX  1 ,•!'•>  1 - 1 ) + ( PU20C  UX1  ,NR1  ) - PH2f  i ( X 1 , ' ■!  K l - 1 ) ) 4 6 R 1 D R 2 
P'f2>>X2  = R'0')(  NX2,  NtJR  - 1 ) + ( PU2P(  NX2  , NR2  ) - PH2i»(r,x2  , ''R2-1 ) )<=DR3PR4 
CPri2iJ=P4  2'.iX  1 +(PH2"K2-P42DX  1 ) *0X10X2  002R°  ) 

I F ( JCHECK  . FD  . <>0  ) i.cmn-CFs  1 t 002900 

PC  I2X  1 = P C i 1 2 ( !<  1 , 1 - t ) + ( °CP2  ( I'  X I , OR  1 ) - PC')2  ( f * 1 , i*m  -1  ) ) 4DR1DR2 

PC'<2  X 2 =PC  )2('.<2,:;P?-1  ) K PC'D?  ( '0  2,  NR  2 ) -PCU2  ( NX  2 , i.  R 2-  1 ) )*Ok3DK4 

CRC-V2  = pri!2  X 1 t-  ( 1 ) *0X1  DX2  0029  3 1 


PU2X  1 = > J2(  I , 1 - n M :J  '12  ( i'!  X 1 , l-  h 1 ) - PVR  (NX  t , UR  1 -1  ) )*n»l  OR  2 
H 1 2 X 2 = R : 2 ( ) < 2 , ' 't  2 -•  1 1 M j J 2 ( Ra  2 , -><■'  2 ) -r  - 2 ( ■*  X 2 , < > 2 - 1 ) ) 4 OH  3 OK  4 


-68- 


CPN2=PN2X1  + (PM2X2-PN2X  1)  * OX  10X2 

P02Xl  = Pn2.(HX  1 ,MR  1 -n  +CP02CNX  1 ,MR  J ) -P02 (NX1 , NR1  - 1)  ) *DR  1 0R2 
P02X2  = PO2(MX2  , LR2-t ) +(PD2(NX2  , UR2  ) -P02  ( N X2.  ,NR2-1  ) )*(>R30R4 
CP02  = P02X1  + (PM2X2-PT2X1  )*0XJI)X2 

PH2Xt=PH2(.MXl  ,‘IRt  ) + (PH2(NXl  ,NR1  )-PH2(NXl  ,NP1-1  ) )*I>H1DP2  _ _ 

PH2X2  = PH2(NX2,MR2-1  ) + ( PH2 ( MX  2 , NR2 ) -PH2 ( NX2 , NR2- 1 ) )*DR30R4 
CPH2=PH2X1+(PH2X2-PH2X1 ) *ox 1DX2 

PCOXtsPCOCNXl ,NRt ) +CPC0CNX1  , NR1-1 ) -PCD (NX l , MR  1 - 1 ) ) *UR 1 0R2  

PCOX2=PCO(HX2,  'OR  2-1  ) + ( PCD  C NX 2 , NR  2 ) -PCO  ( N X 2 , NR  2- 1 ) ) *DR3DR4 
CP.COcPCnXl  + (PCnX2-PCGXl)  *DX1  DX2 

IF(CT.LT.300.  ) C T ~ 3 0 0 . 002940 

IF ( CPH20  . LT . 0 . ) CPH20=1  .G-20 
IFCCPC02.lt. 0. ) CPC02=1.F-20 

IFCCPM2.LT.0.)  C P H 2.  = l , K - 2 0 __  __  __ 

IFCCP02.LT.  0.)  CP02  = 1.£-20 
IFCCPH2.LT. 0.0) CPH2=1 . E-2  0 

I F C C PCO . LT  . 0 . 0 ) CPH2s  1 . E-20_  ___ 

AVERAGE  TEMPERATURES  AMD  PRESSURES  FOR  BLOCK  002980 


IF(J.CT.i)  GO  TO  20  ""  * 

CPH20L=CPH20 
CPC02L=CPC02 
' CPN2L=CPN2 
CP02L=CP02 

CPH2L=CPH2  _ _ ■ __  ' __ 

CPCOL=cpcb “ " * " 

CTL=CT 

20  CQiJTIWUE ' ■ _ _ 

c 

PC02AV=CCPC02+CPC02L)/2.  003000 

PH20AV=CCPH20+CPH2.nL)/2.  002990 

? 2 A V = ( C P U 2 + CP  N 2 L ) * . 5 ’ " . 

PD2AV=CCP02  + CP02M*,5 
PCOAv=ccPcn+cPcor.) /?. . 

PH  2 A V = ( CPH2+CPH2L ) *0.5 ' ~~  


TAV=(CT  + Crij)/2. 

003020 

CPH20L=CPH20 

003030 

CPC02L=CPC02 

CPM2  L=CPM  2 

003040 

CP02  L = CP'J2 
CPH2l,  = CPS2 


1002 

C 


cpcol=cpco  _ 

CTL=CT  ~ 003070 

IF(J.E').l)  GO  TO  1000  003080 

CALL  PAHAMCTAV,  OF.LL  , PH2D  A V , PW2AV  , P 32  A V , PC02 A V ,PCOA  V,PH2AV)  

COOT  I MUR  ‘ 

00  2005  1 1 = t , M W A V ' 003100 

IRLASTC  C n=TRHS(  rn  ‘ 003110 

pc  u )=pc  J ii+iiKcrn 

;)CII)=0(JI)  + AKI|(TT)  0031  70 

IF  CPC.U)  , EO . n.0)  PCTn=1.0E-20  ' ’ 0031  40 

if  cum  .:••).  o.oi  :>cin  = i .of-20  • 003i5j 

A RG=-P(  m / ^OPTr  i . 1 . ?5*P  n T. ) **2/0  c r i ) ) 

IRHSC  l l)=c<P(ARG)  00318'*! 


8 ABC  m = U'>-f,A5T(  rn  - 1'8'ISC  1 C ) } *374  ! 3 . / WAVL5C  l T)  / C FX  P ( 1 . 4 39*’«  A VMO  (1 1 00  3 1 '■*'> 
l ) / T A V ) - l . ) / 7 , 1 4 1 5 R 2 5 1 3 o 003200 

•>  APST8  ci  t ) = ? \ 1 5 .3  l'p  (11)1440(11) 

SRAOC  l 1 ) = -tAi>(  I T)  *1  OOQO.  / -/A  Vi'GC  1 1 ) **2 

SR  AI1.S'/  ( l I )=S3.\r»SVf  ! 1 ) I I ) 

2 0 On  CO  ■>  l’l  I'OE 


l 000  C .1 M r I :ilK  00  3 2 3 0 
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WRITE  (6,202 3)  . - 

2023  FORMAT  ( 20X , ' SPECTRAL  R API ATION ' , / , 1 OX , ' W A VNO ' , 1 OX , ' SR ADSR 1 ) 

00  2022  11=1, MWAV 

2001  sror  = STOT+SR ADSR  (II) 

C wMUTK(6,2024)WAVWO(II)  , SRADSR(IJ) 

c ***TAKE  OUT  WRIT  R STATEMF.NT  OF  SR  ADSR  ( 1 1 ) ***  . __  

2024  FORMAT (2(1  OX, RIO. 4)) 

2022  CONTINUE 

CHECK= l , _ _ 003260 

C 

C RAO( IT)=RADIATTON  FOR  EACH  BLOCK  AT  EACH  WAVELENGTH  (WAVLT(II)  ) 

C . RAO(II)  I N ( W A T T S / C M ? - S R - M T C R 0 N ) ) _ _ _ 

C RAOX(Xl)  IN  ( W AT TS / C M -S R - M T CRON ) ’ 

C TRAOX  IN  (WATTS/CM-SR) 

C RADXT(I)  IN  (WATTS/SR/MICROM) 

C R ADSTR ( 1 1 ) =R ADI  AT  I ON  AT  EACH  X AND  Z SUMMED  OVER  EACH  RATH  LENGTH 


C RADXC I I ) =R ADI ATION  TOTAL  AT  EACH  X-STATIQN  (SUMMED  OVER  Z) 

C TRADX  = RADX(  I)  INTEGRATED  OVER  WAVELENGTH  . _ _003?90 

C RADXTC I)=RADX(I)  INTEGRATED  OVER  ALL  X-STATIONS  ' ~ 003300 

C TOTR AD=THAl)X  INTEGRATED  OVER  X 003310 

C _ _ _ __  „ ' 

1980  CONTINUE 

1990  IF  (M.EQ.l)  GO  TO  1998  003320 

CHECK  = 2 . __ _003  330 

I F ( A L P H A . M E , 0 . . AND. ALPHA .ME. 1 80 . ) GO  TO  1996 

1981  DO  1982  11=1, N WAV  003360 

10=11-1 

i 98  2 RADXT(II)=RADXT(II)  + (RADSTR(II)+RADSTL(II))/2, # PI *DEL 2*4*2* ( ( ZM-1  .) 003370 
l**2-(ZM-2. )**2)  003380 

CHECK  = 3.  . _ _ . _ 003390 

C THIS  CARO  HAS  BEEN  CHANGED  AS  * 

C-  IF  (M.EO.NZPEX)  GO  TO  3001 

IF  (M.EO.NZU)  GO  TO  3001  _ _ _ • 

GO  TO  1998  " 003410 

1996  DO  1997  11  = 1 , MWAV  003420 

C INTEGRATED  RADIATION  IS  DOUBLED  TO  INCLUDE  DOTH  HALVES  OF  PLUME  003430 


RADXC  II  ) =R AOX  (IT)  + (R4DSTR(  ID  + RADSTLC  II ) ) 

C THIS  CARD  HAS  BEEN  CHANGED  AS 

C IF  (M.ED..JZPKX)  RAD(II)=RADX(II)+PADSTR(II)  ... 

IF  (M.EO.NZII)  RAD(II)=RADX(II)+PADSTR(II) 

1997  CONTINUE 

1998  DO  1999  T I = 1 , N W A V _ _ 003450 

C 0=11-1 

1999  RADSTL ( 1 1 ) =R ADSTR ( I [ ) 003460 

2000  CON  T I iNHE  _ 003470 

WRITE  (6,2002) STOT 

2002  FORMAT  (IX, 'SUM  IIP  SR  ADSR  AT  X ST  A TION  1 , 5X  , E 1 0 . 4 ) 

W R I T E ( 6 , 140  2) 

1402  F3RMAT(  8X  , ' NAVE  NUMBER ' , 5 X , 1 W A VF  LENGTH  ' , 6X  , 'RAPT  AT T ON 
$ 2 3 X , ' /i  A V E NUMBER'  ,5X,  'WAVE  LKuGTH  ' ,6X  , 1 RADIATION  ' /) 

WRI TEC6,  14021) 

1 4021  for « AT  ( 1 2 X , ' w / ,s R / C M / M t C R 0 N ' ) 

c 

Ll,NCT  = 0 

DO  12  I 1-1, w WAV, 2 
RADXC  1 1 )=RAI)Xn  I )*f>ri,Z 

f F ( t I . L f . • I A 7 ) RAHXM  T *1  ) = RADX  (T  t+  1 ) *DEL7. 

1F(HADX(  m .I.T.d.  ) v R I I'  E ( 6 , 1 H ) RAOX(Ll) 

18  FIRM  AT  ( ' NEGATIVE  L>  A D l A 1 1 'i  J IS  ILI.IGICIAL  ' ,K14.41 

ir  ( ( 1 1 + 1 ) ..;r . v.jav)  go  n is 

N 3 1 I’E  ( 6 , t „ ) «. .»  v , 1 ( n ) , V A v f,T  ( ] 1 ) , i<  A i > X ( 1 I ) , 

s v a >/<•'.'.(  u n , 'A  vlk  i n- 1 ) , k ft i> x c i r + n 

1 6 E ) R A I (SX  ,F  1 O.D  ,f  1 6.  1 , 7X  , El  2.  •> , ?«X  , F 10 . 0 , F 1 f> . 4 , 7 X , 1 2 . 5 ) 
lit  u: i=li mC r + i 
l F ( I , | .Ci’.I.r.Sl)  GO  T ) 12 
L I 10=0 
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noon  non 


5 


12  CONTINUE 
GO  TO  17 
15  CONTINUE 

WRITE (6, 16)  W A VNO (I I) , WAVLT! II ) ,RADX(II) 
17  CONTINUE 
C 


140  3 FORMAT  ( 1 X , F 1 0 . 0 , F 1 6 , 4 , 7 X , El  2 . 5 ) 0 035  3 0 

DO  2907  I I = 2 , N W A V 003560 

T R A D X = T R A D X + ( W A V f , T ( T I - 1 )-WAVLT(JI)  ) * 

$ (RADXUI-1  )+RADX(II)  )*.5 
2997  CONTINUE 

IF  (I.EO.t)  GO  TO  2999  _ _ 003710 

CENTER  = CEMTER+  (X(T)-X(I-l  ) ) *0 . 5M  TRADX+TR  ADXL  ) *X  ( I ) 

rr)TRAD=TOTRAD+(TRADX  + TRADX(j)*0.5*(X(I)"X(I-l)  ) 00  37  30 

STRAl)(I)  = (TRADX  + TRADXL)*O.S*(X(I)-X(I-l))  _ 

00  2999  1 1 = 1 , N W A V 003740 

RADXTCII)=RAOXT(II)+( ( R ADX ( 1 1 ) +RADXL ( 1 1 ) )/2. 

S )*(X(I)-X(I-t ) )*SC  _ _____  

2999  CONTINUE 

WRITE (6, 50Q0)  TOTRAO 


5000  FORMAT! 1 HO , 24HTOTAL  RADIATION  EM JTTE0= , E10 . 4 , 11H  WATTS/STER)  _ 
WRITE! 6, 1402) 

' DO  5001  1.1  = 1 , MWAV, 2 

IF((.II+1)  .GT.  NW AV)  GO  TO  5002  _ _ _____ 

WRITE! 6,16)  WA VNO ! 1 1 ) , W A VGT ( 1 1 ) , RADXT ! i I ) , WAV N 0(11+1) , WAVLT ( II+i ) 
+ , RADXT! 11  + 1 ) 

5001  CONTINUE  ______  _ _ 

5002  CONTINUE 

WRITE! 6, 16)  W A VNO! IT) , WAVUT! ID , RADXT! II) 


2999  TRAOXMTRAOX 

003760 

DO  2995  1 1 = l , N W A V 

003770 

2995  R A 0 X 0 ! 1 1 ) = R A D X ( 1 1 ) 

093790 

RRN ( I ) =TR ADX 

l 

003790 

WRITE  (6,1492)  I , TR ADX 

OO380O 

1492  FORMAT !/5X,  ’TOTAU  RADIATION  FOR  ST ATION  1 , 1 3 , 1 X , 1 H= , 1 X , El  2 . 5 , 
1 1 WATTS/SR/Cf • ) 


CHECK  FOR  COMBINE  RUM.  IF  SO  OUTPUT  DATA  TO  DISC 

IF  I IiSWTT  . EO  , 2 ) GO  TO  2 2 
WRITE! I DISC)  MWAV, I , T R A D X 

•WRITE!  intso  (WAVK1(  JO)  , >W  AVLT  !v)J  ) ,RADX(JJ)  ,JJ  = 1 ,NwAV) 

T DSCT=  I OSCT  + 1 

22  CONTINUE 

3000  CONTINUE  . _ 003920 

WRITE  (6,1401)  - - 003930 

C E N T R = 0 . 0 003980 

3001  CONTINUE  _ _ _ _ 

3003  CENTRs CENTER  ♦SC/TOTRAD 

TOTRAO=  TOTRAO*,SC  003920  , 

IF  (TOTRAO  .KQ.O.)  T 0 T R A D = 1 ,-*10.  **(-30)  • 

( F ( LSW IT  .FO.l)  WRITE!  I01SC)  TOTR AD , ( R ADX T ( J ) , J = 1 , NW A V j 

3004  WRITE  (6,1404  ) (TTTf,E(J),J=1  ,8  ) 003910 

WRITE  (6,1405) 

1405  FOK'-'Al'  ( 70X  , ’ STATTO'I  RADIATION  (W/SR)1  ) 

WRIT  i‘‘  (6, 1007)  (X(  I n ,PRN(I  I ) .STRAP  (I  I)  , RCAf.C  ( IT  ) , T 1 = 1 , NXS  ) 

1407  FORMAT  (FlS.2,?«X,K13.7,70X,Kl3.?,tX,Al) 

WRITE  f 6,9  121)  004030 

CAINj  F I •)!  A I . j>|,i  IT  mu  T I <1F 


847  1 S'ORW.lT! /I  OX  ,8  M*  DAD!  A ri'l'l  V A 1 1 1 1 F WAS  CAECUM  TED  RASED  ON  THE  SI, OPED  D 1 0 4 0 
4 TvEH.w  THE  P REV!  T I s T *.  0 MOluTS).  004050 

•-RI  TE(R,  Ml  I 1 TOTRAO 

1211  F V 4AT(  l :H> . VOR'li' I'.M.  R A D I A T ! • 1 ■ i Ef.  f T )'  ED  = , E 10.4,11m  * ATTS/STEK  ) 004  OR  0 
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IFCLSWIT  . ED.1)  W R I TE ( I D I SC ) TOTRAD 
W3ITEC6, 1264)  C K M T R 

1254  FORMAT ( 1 HO , 9HCENTROT D= , El  1 . 4 , 1 CM ' ) 

CHECK  FOR  COMBINED  RUN.  SENSE  SWITCH  4 WILL  OVER  RIDE  THE  I, STOP 

t F ( LSwiT .EQ.'l ) GO  TO  24 

GO  TO  (107,103,172) , ISTOP 

DETERMINE  IF  ANOTHER  PASS  IS  TO  BE  PREFORMED 

24  CONTINUE  

JPASS= JPASS+ 1 
IDISCsIDISC+l 

_ KSWITal  ‘ ...  .. 

I F ( J PASS . LT , 2 ) GO  TO  47 
KSWIT=2 

JPASS  = 0 _ _ _ _ 

“ "i  DlSCa  3 ' “**  -----  - --  - - : 

READ  H20  AND  C02  RADIATION  BACK.  FROM  DISC. 

SUM  AND  OUTPUT 

END  FILE  9 _______  __  _ __  _ 

"REWIND  9 ’ 

END  FILE  8 . 

REWIND  8 ___ 

WRITE (6, 27) 

27  F0RMATC1H1 ,30X, ' COMBINED  H20  C02  RADIATI 

$ S U M S ' ) _ _ ■ _ 

W R I T E ( 6 , 2 8 j N W A V 1 , MW A 72  

28  F0RMATC1H  ,/,43X, 'SUMMED  OVER  WAVE  NUMBERS  ’,216) 

DO  2 5 UK  = 1 , NXS 

READ ( 8 ) NWAV,I,CTRAOX 
R E A D ( 9 1 N^AV,I,HTPADX 

CALL  ERRCHK ( 2 ) ' ' 

RE  ADC  8)  (WAVNO(.JO)  , WAVLT(UJ)  ,CRADXCJJ)  ,00  = 1 ,MWAV) 

READ  ( 9 ) (waVMO(J.J)  ,WAVLT(JJ)  ,HRADX(JJ)  , J J=  1 , N W A V ) 

CALL  ERRCHk(3) 

DC)  26  J J = I , N >v  A V 
SRAOX  ( JJ)=HRADX  C.JJ  ) n-CRADX(OJ) 

26  CONTINUE  

HCrOTX ( JK )=HTR ADX+CTRAOX 

N R I T E ( 6 , 9 9 1 1) 

HITEfh  , 1402) 

DO  29  00=1  , UFA V, 2 

IFC  UO+  1)  .GT.UWAV)  GO  TO  30 

WRITE (6, 16)  WAVHO(.U) , WAVLTI JJ) ,SRADX(JJ)  , 

S W A V N f.)  ( J 0 + t ) , w A V L T ( J J + 1 ) ,SRADX(JJ  + 1 ) 

?9  CONI  IN UK 
GO  TO  .11 

30  C INTI" HE 

HlF(h,)h)  WA  VlO(  .1.1)  , 4AVr,T(  JO)  .RRAUXC.JJ) 

31  CONTI  HE 

• U l'  K (6,32)  l , ! I C T O I X ( I < ) 

32  F on  O)  tf // li)v  , » TOTAL  COMBINED  PAD J A I (ON  FUR  STATION  ',T3, 

S ' IS  ' , E l 2 . 4 ) 

K).>'o,o.)  cn  io  vs 

25  c.nrtNNE 

RE 40  (6  ) C I'1 1 < Al)  , ( cn  /,  DV  r ( ) ) , 0=  l , iV.-.'A  V ) 

R F A 0 ( N ) IU'1.14!)  , ( 11  1 O)  / C ( J ) ,J=I  , : v;  A V ) 

CALL  I'  ..  f j 


004090 

OPTION 
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W R I f E ( 6 , 9 9 1 n 
WRITE! 6, 34) 

34  FORMAT  ( 20X , ' W A V N 0 ' , 10X,  ' WAVELENGTH'  , 1QX,  ' HCRADXT  1 //) 

C PUNCH  CARD S FOR  ATMOSPHERIC  TRANSMISSION  PROGRAM 

NNsNWAV 

DO  3 3 00=1  , MWA  V _ 

HCROXT ( 00 ) =HR ADXT ( J I)  +CR  AOXT  ( J j ) 

WRITE(6,35)  WAVNO(JJ) , WAVLT(JJ) , HCROXT (JJ) 

35  F0RMAT(17X,E12.4,6X,E12.4,6XfE12.4)  _ _ 

3 3 C 0 i-l  T I N IJ  E 

S r 0 R A D = H T 0 R A D + C T 0 R A D 

,_WRITE(6f  43)  (TITLE!  J)  , J = 1 , 8)  

43  FORMAT ( t H ///,'  SUMMARIES  FOR  ' , 7 A i 0 ,A2/ ) 

WRITE! 6, 36)  HTORAD 

WRITE(6, 37)  CTORAD  ______  _ 

36  FORMAT  ( 1 OX , ' TOT AD  H20  R ADI ATIOM ' , El  4 . 4 ) 

37  FORMAT! 10X,  'TOTAL  C02  R AD  I ATION  1 , El  4 . 4 ) 

W R I T E ( 6 , 4 1 ) _ _ ___ ___  

41  FORM  AT ! 29X ) . ~ ’ • 

W R I T E ! 6 , 4 2 ) STORAO 

42  FORMAT!'  TOTAL  COM3  INED  R ADI  ATIONJ  , 4X  , El  4 , 4 ) ; ; 

C 

GO  TO  !107 , 107 , 172) , I STOP 
107  CONTINUE 

STOP  ' " ’ ' • 

END  004420 

SUBROUTINE  READER  ___ __  _• . " _ __ 


THIS  ROUTINE  WILL  READ  ONE  COMPLETE  BLOCK  OF 

MO  TOR  DA  T A PER  CAL  L 

COMMON  /INDATA/  X ! 40 ) , R ! 40 , 30 ) , T ! 40 , 30 ) , PC02 ( 40 , 30 ) , 

S PH20(4  0,30) ,PM2I  40,30) ,P02!4O,3O) ,PCO(40,3O)  , PH 2!  40, 30) 
COMMON  /CTROL/  TITLE! 12) , KOU NT , NXS , NORUN 
COMMON  /OTH/  PP , R J , MRP ! 40 ) 

COMMON  /RASP/RASP  _ 

DATA  JJ/1/  “ ‘ ' • ~ 


READ (5,100)  (TtTLE(O)  , J = 1 , 12)  __  __ 

4 R l T E ( 6 , 100)  (TITLE! J) ,J  = 1 , 12) 

100  FJRJMATj  I 2 A 4 1 

A D ( 5 , 101  )NA  Sp)  ^ -CLa.  AcUlcJ-o  y-  . ~ \jc.lku.  w»4- 

READ"!  5,300)  P P , R ■ j , M R P ( j ) ' y 

R J = R 0 M 1 0 Q , 

W«ITE(6,30O)  PP,P,(,M«P(  1) 

300  F3RMAT(2E10.4,I4) 

K01INT  = 0 


NOW  READ  IN  THE  ENTIRE  BLOCK  OF  MOTOR.  DATA 
30  CONTINUE 

READ!  5, 400)  X (JO)  , If.  AST 
■ X ( 00 )=X (00 ) *RJ 

WRITE  ( 6 , 790 ) X ( 00) 

400  FOR..  AT  ( E 1 0 . 3 , 6 9 X , 1 1 ) 

700  FORMAT  ( IX , Et  0. 3) 

!F(  ri.AST.KO.  1)  GO  TO  14 
103  FORM  AT (16) 

H E A 0 ( 6 , 1 03)NtJP(0  I) 

IX--'IKP(OJ) 

DO  12  K K = l , I X 

TKE  PUNCH  DATA  M»UT  FOR  RADI  ATIOM  PROGRAM  |N 
CHANGE  READ  (*J,  !>•»'))  T ( 00  , K K ) , CA'-'D  10 

1 U E E 110.0  «'  i id  ORDER  PPC.J2,  PPtl?d  , PP'j  2 , PPO? 

♦PN2(00,*K)  ,P02M  I.KK) 

READ!  S , SON  ) K(U.I  . x < 1 . f'(  J I . mO  , PC’i'iV  ( J J , KK  ) , I'RVI.If  IJ.KK)  , 
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READ(5,500)R(.L),KK),T(JJ,KK)  ,PCO?( JJ,KK)  , PH2U(,1J,KK) 

500  FORMAT(E10.3,20X,3E10.3)  * j 

RCJJ  ,KK)=R(,)J,  KK)*RJ  I 

PCO(JJ,KK)=l . 0F-20 
PH2CJJ,KK)=1 . 0E-20 

P02(  JO,KK)=t  ,OE-20  

PM2CJJ,KK)=t . 0E-20 

NRIfK(6,f>00)  R(iJ.J,KK),rCJJ,KK)  , PN2C  JJ,KK)  ,PCC1(JJ,KK)  , PC02(JJ,KK), 

+ PH2Q(JJ,KK) rPD2( JJ.KK)  

600  FORMAT  ( IX, RE  10. 4)  ! 

12  CONTINUE  l 

J J = J J + 1 _ ___ ^ 

101  FORMAT  (El  0.4)  ~~  ‘ ~~~  ' ' 

KOUNT— KOUNT+l 

GO  TO  30 ...... _. - 

14  CONTINUE 

CALL,  LI  M I T ( KOUMT , 40 , 1 ) 

CALL  LIMIT (MRP  ( 1 J t 30L2) ' 

RETURN 

_ END  _ __  ....  .J  ' 

SUBROUTINE  l/lMTTC  I ARG,LIM  ,NPS) 

c * 

IFCIARG.GE.O. AMD  . I ARG . LE . LI M ) RETURN __  ' 

c " ” - . - ■ - • 

4 R I T E ( 6 > 1 0 ) I A R G , L I M , N P S 

10  FORMAT  C 1 ILLEGAL  INPUT  PARAMETER  • , • ! 

$ 2110/'  CALLED  FROM  LOCATION  ' , I R ) 

C 

___  STOP ; 

END  • 

SUBROUTINE  SWAP 

C 

COMMON  /BAN  OP/  XKK(254,7)  , DD.T  NV  ( 254 , 7 ) , UK ( 254 ) , AKU (254 } 

COMMON  / K P M 1 / N W A V , W A V M 0 ( 2 5 4 ) , W A V L T ( 2 5 4 ) 

C _ 

DIMENSION  X SAVE (254, 7) , OSAVEC 254 , 7 ) 

C 

DO  10  00  = 1,  NWAV 
DO  10  K<  = 1 , 7 
X H 0 L 0 = X K K ( 00  , K K ) 

D H 0 L 0 = 0 D l N V ( 0 0 , K K ) _ ..  . 

X KK  ( 0 J , K K ) = X S A VK  ( J J , KK  ) 

D !)  I N V ( . ! 0 , K.  K ) r o S A V E ( J J , K K ) 

XSAVE(OJ,KK)=XHOLD 

DSAVE(.IJ,KK)=DHO[,D 

10  CONTINUE 

c 

NWVOLI)=NWAV 
si  .‘1 A V =.  N 'M,ST 
NWLST  = U 4V0L!) 

C 

RETURN 

C 

ENTRY  SAVE1 
C 

DO  11  ,hl=),NWAV 

)0  11  K*=l,7 

X S A V E ( . 1 , 1 , K K 1 = X K K ( J J , K K ) 

■ ) S A V E ( . I J , <K)=r.i)r  Vtf(  J ! , <<  ) 

11  CO  "j  Tl'IilK 
C 

)VLST=NRAV 

c 

ne  r s i »•>  i 

Emu 
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SUBROUTI OK  PRMD 

COMMON  /INOATA/  X ( 4 0 ) , R ( A 0 , 30 ) , T ( 4 0 , 3 0 ) , PC02 ( 4 0 , 30 ) , 

S P H 2 0 ( 4 0 , 30) , PH? (40, 30) ,PH2(40, 30) , PCO(40, 30) ,PH2(40,30) 

COMMON  /CTHOL/  TITLE ( 1?)  , KOUNT , NXS , NORUN 
COMMON  /OTH/  PP,.RJ,NRP(  40)  

WRITE(6,10) 

.10  FORMAT!  1 HI ) 

WRITE (6,11)  NOR UN 

11  FORMAT  ( 1 H , 2 5 X , 1 MOTOR  *,13//) 

__  WRI  TE(6, 12).  KQtlMT  

12  FORMAT!  I H'  , 9X , 1 NUMBER  OF  X STATIONS  IS  1 , 14) 

DO  14  .JiJ  = l , KOUMT ___  _ 

NRP1=NRP( JJ) 

W R I T E ! 6 , 15)  JJ,X(JJ) ,MRP1 
J_5_.F?«>lA:r(//.'  X ( ' , t 2 , 1 ) = 1 , E l 2 . 4 , 

$ 8X  , ' NUMBER  OF  POINTS  AT  THIS  X STATION  IS'  ,14) 

WRITEC6, 16) 

16  FORMAT!  1H  , 4 X , ' R ' , 3 3 X , ' T ' , 11  X , ' C02  * , 1 OX , _ _ 

$ ' H 2 0 ' » 1 1 X , 1 N 2 1 ,1 IX, '02'/) 

00  14  KK=1,NRP1 

.WRITE! 6, 18)  R(JJ,KK) ,T( JJ,KK) ,PC02(JJ,KK)  , 

S PH20 ! J J , KK ) ,PN2(JJ,KK) ,P02(JJ,KK) 

18  FORM  AT ! 1 H , 6 ( E 1 2 . 4 , 1 X ) ) 

JJ  CONTINUE 

RETURN 

END  __  _ 

subroutine  bpprt 

THIS  ROUTINE  PRINTS  THE  9AMO  PARAMETERS  TOGETHER  

■WITH  THE  ASSOCIATED  WAVENUMBERS  AND  WAVELENGTHS 

COMMON  / B A N 0 P / XKK ( 2 54 , 7 ) , DPI N V ( 2 64 , 7 ) , UK ( 254 ) , AKU ! 254 ) 
COMMON  / R P W 1 / N W A V , W A V N 0 ! 2 5 4 ) , W A V L T ! 2 S 4 ) 

C 

NRITE!6,U) 

11  FORMAT!  1H1  ,S1X,  "3  A At  0 P A R M F.  T E R S'//) 

L1NCT=0 
WRITE (6, 12) 

t 2 FORMAT! 1 OX  , ' WAVENUMBER ' , 1 X , ' WAVELENGTH  ' , 20 X, 

S 'XKK  AND  OOINV  PARAMETERS'/) 

C 

DO  10  iM  = l , MWAV 

WRITE! 6, 14 ) WAV NO!  J J ) , W A VLT ! JJ ) , (XKK!  JJ,K) ,K  = l ,7) 

WR  ITE(6 , 14  ) WAVNO!  JO)  , *‘AVLT(vU  ) , ( Df)  f W V ! U J , K ) ,-Kr.l  , 7 ) 

14  FORMAT! 14X,FS.0,4X,KR.l , 7X,7E12.4) 

WRITE (6, 16) 

15  FORMAT!  lit  ) ' 

LINCTsUNCT+3 

IF! LIMCT.LT.561  GO  TO  10 
LI JCT=0 
WRITE! 6, 16) 
lb  F OR M A T ! I H 1 ) 

WRITE (6, 12) 
to  comuuK 
c 

R TURw 
EvO 

SUM  rod  r ( JK  r.  f 1/  (TS(  ARRAV  , NO,  M.Mfc  , M<pPEK  ) 

C 

R <*.•  A 1 1 A M A V (254) 

C 

R H l 0 M = 1 . E “ <0 
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DO  10  J,J=l,NO 

IF  ( ARRAY  (JO)  . GT.RHIGH)  RHIGH=ARRAY  ( J J ) 
10  CONTINUE 
C 

IUP  = A1.0G10(RHIGH)  + 1 

ILP=IUP-3  

RUPPER=io.**TUP 

RL0W=10.**n,P 


RETURN 

END 

__  SUBROUTINE  ERRCUKCTl.OC) _ 

RETURN 

END 

SUBROUTINE  PAR  AM  ( TAV  , X L , PH.20  , PN2  , P02  , PC02  , PCO , PH2)  

COM MON /CHECK/ 1 READ 

COMMON  /BAMDP/  XKKC254,7) , DDINV ( 254 , 7 ) ,UK(254)  ,AKU(254) 

COMMON  / R P N t / MWAV,NAVMO(?54)  , W A VLT  (.2  54.) 

“common  /snitch/  kswir  , ixxjg 
real  T T ( 7 ) 

REAL  LOG  l , LOG2  _ 

C 

DATA  TT(1  ) , TT  ( 2 ) , TT ( 3 ) ,T1’  (4)  , TT ( 5 ) , TT ( 6 ) ,TT(7)/300.  ,600.  , 1000,  , 

SI  500.  ,2000,  ,2500.,  3000,/ 

DATA  I0BG/-1/  004470 

T=TA  V 

IXXJGsIXXJGt  1 ..  ..  ; ' ' ' 

IF  ( K S W IT. WE. 2)  GO  TO  23 
T T ( 3 ) = 1 2 0 0 . 

TT(4)  = 1SOO. . 

TT( 5 ) =1 800 . 

TT ( 6 ) =2400 . 

23  CONTINUE  _ _ __ 

1DBG=I0»G+1  ■ 004480 

C CALCULATE  OPTICAL  DEPTH,  (J,  AND  WIDTH,  GAMMA  004490 

C _ 

1273=273. /T 
S0T27  3=S0RT  f T273 ) 

I F ( K S W I T , E 0 . 2 ) GO  TO  18  _ ; 

U=T273*PH20*XL 

GAMMA  = t 20. 1 *PH?.0/TtS3T27  3*(  . 09* ( PH20+PN2 ) + . 0 4*P02+ . 1 2*PCI32 
$ +.05*PH2+0.1*dC0) 

GO  TO  19  " 

18  CONTINUE 

U =T  2 7 3 * PCO  2 * X L 

GAMMA  = 0.0J  *PC02*T273+(  . 07  *PN  2 + 0 . 055*P02+0 . 09*PC02+PH2*  . 08 
S + PCO  * 0 . 0 6 ) * S 0 T 2 7 3 

19  CONTINUE 


SELECTION  OF  TEMPERATURES  TO  6E  USED  IN  INTERPOLATIONS  004530 

t = l 

DO  30  KK  = 1 ,7  

K = KK 

TT1  = AHS( T-TTf K ) ) 

IF  ( TTI  .I.E.E-20)  GO  TO  31 

30  CONTINUE 
GO  ro  32 

31  = 

111  17  = 001  iVn  , K 1 
GO  M 16 

32  com  •NIK 

n 0 = 3,7 

X = I 

1 r ( f I’  C 0 ) ” M R , K , 9 
3 CHI' i M.E 
I C I I f I 'll,.; 
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IFC T.GK,TT(K-1 ) ) K = K + 1 
TTK=TT(K)-TT(K-2) 

TTKlsTT(K) -TT(K-l  ) 

TTK2=TT(K-l  )-TT(K-2) 

DI V=TTK*TTK  1 *TTK2 

C ' _____ 

r'  1000  '"CONTINUE ~ 

LDG1  = ALOG (XKK(I,X-1  )/XKK(l,K~2)) 

l.OGX  = Al.nG(XKKCt,K)/XKK(  t,K-2))  

6=(  (TTK**2)*L0Gt  -(TTK2**2)*LnG2)/DIV 
C- ( TTK2*  L0G2 -TTK*L  OG l ) /OIV 

XK=EXP( ALOG (XKK ( T ,K-2})+R*(T“TT(K-2))+C*(.r’*TT(K-2))**2_)' 

C 

IF(T  .IjT ,TT(2) ) GO  TO  14 
0 LOG = A LOG (DDIMV(J,K-2)) 

DINV=EXP(  ( ALOGCnOPIVC  I , K-l  ) ) -DLOG ) * (T-TT  ( K-2 ) ) /(  T T ( K - 1 )-TT(K-?.)  ) 
C +DLOG) 

GO  TO  16 
C 

14  CONTINUE 

OLQG  = ALOG ( DD  IN V ( 1 , 1 ) ) 

l)  I N V = EXP  ( ( ALOG  (ODIN  V ( 1 ,2  ) ) -DLOG)  * ( T-TT  ( 1)7  / ( TT  ( 2 ) -TT  ( 1 ) ) + D L 0 G ) 

16  CONTINUE 


CALCULATE  PARAMETERS 
UK(I)=XKtIJ 

AKU(I)=GAMMA*DINVti.IK(I) 

IF  C T T 1 . L E.  E - 2 0 1 G Q_  TO  1 7 

KDEX=K-2 

ioex=k-i 

FL0'N=:.9*AMIM1 (XKK(T,KOEX) ,XKK(I, IDEX) ) 

FHIGH=1 . l*AMAXl (XKK(T,KDEX) , XKK (I , IDEX) ) 

IFIXK.GE. FLOW. AMO. XK.LE.FHIGH)  GO  TO  17 
'6  = ALOGCXKKCl,K-ll/XKK(T,K-2. 11/TTK2 
XK=EXP ( ALOG  ( XKK  ( T , K-2 ) ) +B* (T-TT ( K-2 ) ) ) 

IFCXK.LT. FLOW . OR.XK.GT. E H I G H ) W RITE (6,20) 

$ WAVUO(I)  , I , XK  , T , II , GA  MM  A , T'i'K. , 0 , C , ()  IN  V , UK  ( I ) , AKU(I) 

20  FORMAT (IX,/, 1 XK  STILL  OUT  OE  RANGE  AT  WAVNO  = ',E10.0,'  AMD  I = 

S ,I5,/,1X,5EJ4.4/,1X,5514.4) _ __ 

17  CONTINUE 
1 = 1+1 

IF  ( I-N  /IAV)  29 , 29, 1 5 __  __ 

29  IF  CTT1  .LE.K-20)GO  TO  31  “ ” ~ 

IFCI-tMAV)  1000,1000,15 
C 

1.5  CONTINUE  ” 

I RKAI)=  I READ+ 1 
RETURN 

END  

/ / G 0 . F r 0 3 F 0 0 1 01)  OSN  = HARNE.PLOTPf uni t=sysd a, disp=( new,  delete ) , 

//  OC-3=(RECFM  = VS«,LRECL=1  00,8LKSTZE=1000,6UFNO=1)  , 

//  SPACE=(CXL,  (5,5)  , I’LSE) 

//GO  . F TOOK  001  0 0 OSNsMAR  •'  'f  . P 1,0  TO  , U N.IT  = S Y SO  A , OISP=(  WF.W  , OF  LFTE  ) , 

//  IJCi>=(PKCF  ''  = VSi<  , LRECI.  = 1 00, 8LKS  IZE=  1 000  , fU)FoO  = l ) , 

//  SP aCK=(C CL, (5,5) , RLSF ) 

//GO. SI  SIN  OH  * 
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INPUT  INSTRUCTIONS  AND  SAMPLE  INPUT  DATA 


APPENDIX  II 


The  computer  program  input  is  composed  of  two  parts:  the  flow  field  section 
and  the  band  parameter  section.  The  flow  field  of  a plume  can  be  generated  by 
several  computer  codes  like  the  NASA  MOC  code  or  the  AeroChem  code  [3  and  5]. 

The  flow  field  program  provides  the  temperature  and  partial  pressure  at  different 
locations  along  the  plume.  This  information  is  necessary  to  calculate  the  optical 
depth  and  the  radiance  of  the  plume.  The  band  parameter  input  is  composed  of  the 
coefficient  of  absorption  part  and  the  line  density  part.  It  is  tabulated 
according  to  the  wave  number  and  temperature.  Subroutine  READER  reads  in  the 
flow  field  data  and  subroutine  PARAM  calculates  the  transmissivity  T. 

Table  A.l  lists  all  input  parameters,  format  of  the  read  statements,  and 
their  meaning,  in  the  same  order  as  encountered  in  the  program  for  the  flow  field. 
Table  A. 2 lists  the  input  parameters  required  for  the  radiation  calculations. 
Tables  A.l  and  A. 2 indicate  all  input  data  needed  for  operation  of  the  code. 

A listing  of  a sample  input  data  case  is  included  as  Table  A. 3. 
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TABLE  A.  1 


Card  # 

Parameter 

Format 

Meaning 

1 . 

Title 

12A4 

Title  of  the  Flowfield 

2 

PP 

E10.4 

Total  Pressure 

RJ 

E10.4 

Radius  of  the  Nozzle 

NRP(l) 

14 

//  of  Radial  Points  at  the 
Beginning  of  the  Plume 

3 

X(JJ) 

E10.3 

Axial  Location 

ILAST 

69X.I1 

Control  value.  ILAST  = 1 
for  the  end  of  the  flow 
field 

4 

NRP(JJ) 

15 

# of  Radial  Points  at 

Axial  Location  X(JJ) 

5 

R(JJ,KK) 

E10.3 

Radial  Location  at  Point 
(JJ,KK) 

T(JJ,KK) 

E10.3 

Temperature  at  Point 
(JJ,KK) 

pco2(jj,kk) 

E10.3 

Partial  Pressure  of  C07  at 
Point  (JJ,KK) 

ph2o(jj,kk) 

E10.3 

Partial  Pressure  of  I^O  at 
Point  (JJ.KK) 

Repeat  card  5 for  NRP(JJ)  times. 

Go  back  to  card  3 for  another  axial  point. 

At  the  end  of  the  flow  field,  set  HAST  = 1.  Such  value  of  ILAST  will 
terminate  the  reader  subroutine  and  return  to  the  MAIN  ROUTINE. 


TABLE  A.  2 


Card  # 

Parameter 

Format 

Meaning 

1 

NWAV  1 

110 

Wave  number,  lower  limit 

NWAV  2 

. 110 

Wave  number,  upper  limit 

2 

WAV  NO (I) 

F9.1 

Wave  number 

XKK(I, J) 

7E9.2 

Coefficient  of  absorption 

K <§  300,600,1200,1500, 
1800,2400,3000  °K  of 
wave  number  WAVNO(I). 

3 

WAVNO(I) 

F9.1 

Wave  number  (must  be  the 
same  as  WAVNO(I)  in  card  2) 

DDINV (I, J) 

7E9.2 

Line  density  at  temperature 
of  300,600,1200,1500,1800, 
2400,3000  °K  of  wave 
number  WAVNO(I). 

Repeat  card 

2 and  3 for  N = (NWAV  2 

- NWAV  l)/25  +1 

4 ALPHA  F10. 5 Aspect  angle  of  the  line  of  sight 

ISTOP  15  ISTOP  = 1,  will  terminate  the  program 

after  single  alpha  calculation.  ISTOP  = 2 
will  return  to  this  point  for  another 
ALPHA  calculation.  ISTOP  = 3,  will  read 
in  a new  set  of  flow  field. 
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TABLE  A. 3 


.10006  010.39306  00  1? 

0.0 


0.0  0.24006 

0. 39306-010.24006 
0.78606-010.24006 
0.11796  0 0 0 . 2 4 0 0 £ 
0.15726  000.24006 
0.19656  0 00.2400 £ 
0.23586  000.24006 
0.27516  000.24006 
0.31446  000.24006 
0.3537E  00  0. 2 4 00  E 
0.39306  000.24006 
0.4606E  000. 2880 E 
0.4106  00 
0.0  0.24776 

0.39936-010.24776 
0 . 7 9 8 2 E - 0 1 0 . 2. 4 7 7 E 
0.11976  000. 247 7 E 
0.15966  000.2477 £ 
0. 19966  000. 247 7 K 
0.23956  000.24776 
0,27945  000.24776 
0 . 3 1 93r  000.24756 
0.3591c.  000.24486 
0.39976  000.22286 
0.45316  000.12126 
0.54396  000.35206 
0.64496  000.35396 
0.74596'  000.30206 
0.83636  000. 28806 
0.8146  00 
0.0  0.24876 

0. 39976-010. 24876 
0.79946-010.24876 
0.11996  000.24876 
0.13995  000.248/6 
0.19996  000.248/6 
0.23986  00 0.24876 
0.27986  000.24866 
0.31976  000.24736 
0.35986  000.2.4126 
0.40166  000.21736 
0.45396  000.14/96 
0.533/6  000.81666 
0 , o 3 0 8 6 0 0 0.30 i 2 6 
0.73186  000.37396 
0.82696  0 00.32316 
0.91886  Oo o.J 017 6 
0.1  0036  0 1 0 . 7 v 2 3 6 
0.1226  01 


0 4 0. 1000 6 -1  90. 10006-1 9 0 . 75466-010.25636  000.10006-19 
040 . lOOOE-1 90 . 1 0006-190. 75466-010.25636  000.10006-19 
040 . 10006- 1 90. 1 0006-190. 75466-010.25636  000.10006-19 
040. 10006-1 90. 10006-190. 75466-010.25636  000.10006-19 
040. 100 OE -190. 10006-190. 7 5466-010. 2 5636  000.10006-19 

040.10005- 190. 1000E-190 . 75466-01 0. 25636  000.10006-19 

040. 10006- 1 90 . 1 0006- 190. 7 5 466-0 10. 25536  000.10006-19 
0 4 0.1  000 6-1  90. 10 OOF.-  1 9 0.75466-010.25636  000.10006-19 

040.10006- 190.10006-190.75466-010.25636  000.10006-19 

040. 1 0006- 190. 10006-190. 75466-010.25636  000.10006-19 
040.10  00  6-190.10006-190.73466-010.25636  000. 1000 K -19 

030.10006- 190.10006-190.28836-060.28656-060.10006-19 

040. 10006- 190. 10006-19 0.7 9546-0 10. 2625E  000.10006-19 

040.10006- 190.10006-190.79546-010.26256  000.10006-19 

040.10006- 190. 10006-190. 79546-010.26256  000.10006-19 

040. 1 0006-  l 90. 10 006" 190. 7 9546-0 10. 2 6 256  000.10006-19 
0 4 0 . 1 0006  - 1 9 0 . 1 0 00  F.- 1 9 0 . 7 9 5 4 b - 0 1 0 . 26  256  000.1000  6-19 
0 40. 100 06-1  90. 10006-190. 795 46-010. 26 2 5 6 000.1000  6-19 

040. 10006- 190. 1 0006- 190. 7 9346-0 1 0 . 26256  000.10006-19 

0 4 0. 1 000 6-1  90. 1000E-1  9 0.7953c. -010. 26256  000.10006-19 

040.10006- 190.10006-190.79476-010.26236  000. 10006-19 

040. 10006- 190. 10 006-190. 7 83 46-010. 2592b  000.10006-19 
040.  10006.-1  90. 100UF.-1  90. 6 8446-01  0.22796  00  0.10  006-19 
040. 1000E-190. 10  OOF- 190. 2 7 096-010. 84656-010, 10006-19 

030. 10006-  190. 100 06-1 90.0866 6 -020. 2061 6 -010. 10006- 19 

030 . 1 0006-1  90. 1 0006-190 . 16606-020. 4 98b£-020 . 100  06-19 

030. 10006-  190. 10 006- 190. 34o76-0 3 0.1 0676-020. 10006-1 9 

030.10006- 1 90. 10006- 190. 288 36-0 60. 288 56-060. 1000E- 19 

0 40.  10006-1  90. 1000F- 1 90, 82096-0  10,260951  000.100  06-1  9 
040. 100 OK- 1 90. 10006-190. 8 2 09b- 01 0.26096  000.10006-19 
0 40. 10006-190. 1000 6- 190. 82096 -010. 2609  6*  000.10006-19 
04  0 . 1 0005-  1 90. 1 0006-1 90. 8 2 096-010 .26096  00 u . 10006-19 
040. ] 0006- 1 90. 10006-190 . 82096-01 0.26096  000.10006-19 

040. 1 noot’-l  9 0 . 10006-1  90 . 82  0 96-01  0.26096  000.10006-19 
0 4 0 . 1 0 0 0 6 - 1 9 0 . 1 0 0 0 r.  - 1 9 0 . 8 2 0 9 6 - 0 1 0 . 2 6 0 9 b 0 00.1  0006-19 

040. 10006- 1  90 . 1000b- 190. 8 2 07 b-0  10.260  7t.  00  0.10006-19 

040. 10006- 1 90.10006-190.81856-010.25946  000.10006-19 
04  0. 100 06-1  90. 10  006-1  90. 80  356-010 .2512b  000.10006-19 
0 40. 1 0006-1  90  . 10006-1  90. 7 2 12.6-01  0.21  3 7b  000.100  06-1  9 
o 4 0 . 1 0 0 0 6 - 1 9 0 . 1 0 0 0 6 - 1 9 0 . 4 2 6 1 6 - 0 1 0 . 1 2 2 0 1 0 0 0.10006-19 
o 3 0 . 1 0 0 0 £ - 1 9 0 . 10006-  190.1  3 8 lb-010. 4558  6 -01 0.1000  6-19 
0 30. I 00  06-1  90 . 1 0006-1 90 . 57 4 7 b- 02  0 . 1 6876-0 J 0 . 1 0006- 1 9 
0 30. 1040  6 -1  9 0.  10006-1 9 0.2  2 40  F -020. 83  7 16-02  0.100  06-19 
0 40. 10006-1  90. 100 OF- 1 9 0.90  2 96-03  0 ,270  8c -020 . 1000b- 19 
>'  10.  I (.  006-1  70. 10 OOF-1  90. 3 301  b-03  0 . 105  36-020. 100  06-1  9 
0 30. 1 0006-  1 90. 1 000c.- 1 90 . 1 1 48  6-0  3 0 . 34D4b-030 . 1 00  06-  1 9 


0.0  0.249 7 K 

0. 40055-010. 2-.V/C 
0 . 8 0 U 9 b - 0 t u . 2 i v / F 
0.1201c.  o o o . 2 1 9 7 c. 
O.loOlb  000.24476 

0.2  )0 2b  0 0 0 , 2 t 2 / 6 

0.24026  00 o.2443c 


0 4 0 . t GOOF-t  90.  1 0 00b- 1 9 0 . 8 4 020010.259  86  000  . luOOE-l  9 

0 40  . t 0006-  1 90 . 1 0006- 140 . 84026-0  1 o . 23986  000 . 1 00  36  - 1 9 
04 0. 10006-1 90. 1 0006-190.84026-01 0.2598C  000. lOOOb-) 9 
•'4  0. 1 on  >6  - 1 40. 1 0 0 06-  1 9 0,8  4 026-010 .23986  000  . 100  Ob- 19 
<’40.  I 0006- 1 90 , 1 OOUb- J 90.84026-0 1 0.25486  000.1000f.-19 

1 4 u , 1 OO.jF-1  9 0 . 1 UO  U C - 1 9 0 . 8 « 0 2 6 - 0 1 0 . 2 5 4 8 b 0 0 o . 1 0 0 0 6 - 1 9 
>'4  0 . 1 0 0 0 6- 1 9 0 , 1 0 U 0 6 - l 9 0 . 8 4 0 4 r, - 0 1 u .? 3 9 o 6 oOO.  1000b- 19 


> 
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0 . 2 8 0 2 E 0Q0.24H9E  040 . 1 OOOE-  HO  . I OOOE-  1 90  . H 4 1 1 K-0  1 0 . 25m  8E  000.1000E-19 
0.3202b  000.24b3E  040 . 1 00 OE- 1 90 . 1 OOOE- 1 90 . H 4 20E-0 1 0 . 2550E  000.1000E-19 
0.3608E  0 0 0 . 2 3 a b E 040  . 1 00  0E  - 1 90 . 1 000t>  1 90  . d 3 46E-0 1 0 . 2 4 1 3E  000.1000E-19  ' 

0.4052E  000. 2 161 E 040 . 1 000E- 1 90 . 1 OOOE- 1 90 . 7 7 1 8F>0 1 0 . 20 32E  000.1000E-19 
0.4602E  000. 1643  E 040 . 1 OOOE-  1 90 . 1 000b-  1 90 . 5494E-0  l 0 . 1 359E  000 . 1 00  O b- 1 9 
0.5349b  000.1050b  040.1000b-190.1000E-190.26bOh-010.b755E-OlO.1000fc.-19 
0.62  53E  OOO.b'/OOE  030 . 1 OOOE-  1 90 . 1 OOOE-  1 9 0 . 1 1 b4E-0 1 0 . 30b9E-0  1 0 . 1 OOOE-  1 9 
0.7239E  000. 478 JE  030 . 1 00 OE- 1 90 . I OOOE- 1 90 . 63 60 K-0 20 . 1 466E-0 1 0 , 1 OOOE- 1 9 
0.8210E  000, 3343 K 030 . 1 00  Ob  - 1 90 . 1 OOOE- 1 90 . 26 30E-0 20 . 7 40 OE-020 . 1 00 OE- 1 9 
0. 91 45E  000. 3379 E 030 . 1 00 OE- 1 90 . 1 UOOE- 1 90 . 1 3 39E-0 20 . 3H 4 2E-020 . 1 OOOE- 1 9 
0.1004E  010. 313b £ 0 30 . 1 00 OE- 1 90 . 1 OOOE- 1 90 . 6796E-0 30  . i 97 4E-0 20 . 1 0 OOE- 1 9 
0.1009E  010. 3 004 E 0 30. 1000 £-190. 1000E -190. 32 6 5E-030. 9 5o7E-030. 1000 E- 19 
0. U72E  010. 29 29 E 0 30 . 1 00 0 E- 1 9 0 . 1 OOOE- 1 90 . 1 2 86 E- 0 3 0 . 37 8 bE-03 0 . 1 00 OE- 1 9 
0.161E  01 

0.0  0.2506E  040. 1 OOOE-1 90 . 1 OOOE- 1 90 . 8 543E-0 1 0 .259 1 E 000.1000E-19 

0 . 1 0 0 9 E - 0 1 0 . 2 3 0 b E 040.1000E-190.1000E-190.8543E-010.2591E  000.1000E-19 
0.8017E-010.250ob  040 . 1 00  OK  - l 90 . 1 OOOE-1 90 . 8 543 E-0 1 0 . 259 1 E 000.1 OOOE- 19 
0 . 1 2 0 3 £ 000.2506E  040 . 1 OOOE- 1 90 . 1 OOOE-1 90 . 8543E-0 1 0 . 259 1 E 000 . 1 000E-19 
0.1604E  000. 2503  E.  040 . 1 00  OE- 1 90 . 1 OOOE- 1 90 . 8 544E-0 1 0 . 259  1 E 000.1000E-19  * 

0.2004E  000. 250  4E  040 . 1 OOOE- 1 90 . 1 OOOE- 1 90 . 8 549E-0  l 0 . 2589k.  000.1000E-19 
0.2405E  000.2500b'  040 . 1 OOOE- 1 90 . 1 OOOE-  1 90 . 85  ou£-0  1 0 . 2583b  0 0 0 „ 1 0 0 OE- 1 9 
0.2a05£  000. 2488  E 040 . 1 OOOE- 1 90 . 1 OOOE- 1 90 . 86U9E-0 10 . 2562E  000.  l'OOOE- 19 
0.32U8E  000.2456E  040 . 1 OOOE- 1 90 . 1 OOOE- 1 90 . 8 o87£-01 0 . 2496E  000.1000E-19 
0.3625E  000. 23 7 5£  040 . 1 OOOE- 1 90 . 1 OOOE- 1 90 . 8 70 1 E-0 1 0 . 23 25E  000.1000E-19 
0.4087E  000.2168E  0 40 . 1 OOOE- 1 90 . 1 OOOE- 1 90 . 8206E-0 1 0 . 1 964E  000.1000E-19 
0.4649b  000. 1746 E 0 40 . 1 0 0 Ob- 1 90 . 1 OOOE- 1 9 0 . 64 01 b-0 1 0 . 1 4 1 6£  000.1000b-l9 
0.5368E  000. 122b£  040 . 1 OOOE- 1 90 . 1 OOOE-l 90 . 3708E-01 0 . 834b£-0 1 0 . 1 OOOE-1 9 
0.6237E  000. 826  7 E 030 . 1 OOOE- 1 90 . 1 OOOE-l 90 . 1 8 48E-0  l 0 . 4325E-0 1 0 . 1 0 00b- 1 9 
0.7196b  000.5889b  030 . 1 OOOE- 1 90 . 1 OOOE- 1 90 . 934ob-020 . 228 4E-0 1 0 . 1 OOOE- 1 9 
0.8170b  000. 45 90 E 030. I OOOE- 1 90 . 1000E- 1 90 . 5004E-020 . 1 27 1 £-01 0 . 1 000b- 1 9 
0.9122E  000.3867b  0 30 . 1 OOOE- 1 90 . 1 OOOE-1 90 . 2a 1 3 E-0 20 .7 375E-02 0 . 1 OOOE- 1 9 
0.1003E  010.3461E  0 30 . 1 OOOE- 1 90 . t OOOE- 1 90 . 1 62 7E-020 . 4369E-02 0 . 1 000b- 1 9 . 

0. 109 lb  010. 32  20t  030 . I OOOE- 1 90 , 1 OOOE-l 90 . 94 1 4E-030 . 2 57 2E-02 0 . 1 OOOE- 1 9 
0.1175E  010. 30 7 2E  0 30 . 1 OOOE- 1 90 . 1 OOOE-1 90 . 5280E-030 . 1 4b 2 k. -02 0 . 1 OOOE- 1 9 
0.1  25  7E  01  0.298  OR  030 . 1 0 0 OE- 1 90 . 1 OOOE-l  9 0 . 2.741  E-0  30 . 7655E-030 . 1 OOOE-  1 9 
0.1337b  0 1 0 , 2 9 2 2 £ 030.1000E-190.1000E-190.1138E-030.3194E-030.1000E-19 
0.197E  01 

0.0  0.2513b  040. 1000K-1 90. 1 OOOE- 1 9 0 . 8 6 39E-0 1 0 . 258 oE  000.1000E-19 

0. 8029E-O1 0.251 3£  o 40 . 1 OOOE- 1 90 . 1 OOOE- 1 90 . 8 640E-0 1 0 . 25aaE  000. 1000b- 19 
0.1603E  o 0 0 . 2 3 1 2 E 0 40 . 1 000 E - 1 90 . 1 OuOk.- I 90  . a 6 4 7 E-0  1 0 . 258 bb  OoO.lOOOE-19 
0.2407  b 000. 230  3 k.  0 40 . 1 0 0 0k.- 1 90 . 1 OOOE-l  90 . 8 7 0 8 E-0 1 0 . 23b  8 b 00  0 . 1 OOOE-1  9 
0.3218E  000.2  154b  0 40  . ].  OOC  E-  1 90 . 1 OOOE-  1 90  . d 9 37E-0 1 0 . 2444 1 000.1000E-19 
0.4126b  0 0 0 . 2 1 7 o E 1)40 . 1 OOOE- 1 90 . 1 OOOE-1 90 . 855  aK-01  0 . 1 91  ob  000.1000E-19 
0.5388b  OUO.  15  56E  040 . 1 OOOE-  1 90 . 1 OuOE- 1 90 . 4362  E-0 1 0 . 93H7E-0  .1  0 . 1 OOOE-  1 9 
0.7168E  000.69  35c,  030 . 1 00  OE- l 90 . 1 OOOE- 1 90 . 1 57  7 E-0  1 0 . 3044E-0 1 0 . 1 OOOE-  1 9 
0. 9102b  000 , 4 412 E 0 3 0 . 1 OOOE- 1 90 . 1 uOOfc - 1 9 0 . 4 6 3 1 E-0 2 0 . 1 1 o 4E-U 1 0 . 1 OOOE- 1 9 
0. 109  lb  0 10.3494b  030. 1 oOOE-1 90. 1 OOOE- 190 . 1 787E-020.4492fc.-020. 1 OOOE-1 9 
0. 12oOE  010.  3120c,  0 30. 1 OOOE-1  90. 1 OOuE- 1 90 . 6 8 57  E-0  30 . 1 7 8 1 E-0  2 0 . 1 OOOE- 1 9 
0.1120E  01 0.294  IE  030. 1 000 E-190. luOO E -190.1 7 32E-030. 4 571 E-030.1000E-19 
0.237E  01 

0.0  0.2520  r,  0 4 0 . 1 9 0 3 £ - 1 9 0 . 1 0 0 0 k.  - 1 9 0 . 8 7 2 1 E - 0 1 0 . 2 5 6 5 E 000.1000E-19 

0 . 8037c.-  J 10. 7524k.  0 40  . i i,00t-  1 90 . 1 OOuE- j 9 0 . 87  25b-0 1 0 . 2585b  000 . 1 OOOE-  19 
0.  1607  k.  005.26l.-E,  0-*0 . 1 KOOb-  1 90 . 1 OOOE-  1 90  . b/bJE-OlO.  2576b  00O.1OO0E-19 
0.241  1 E j o 0 , / 5 " 6 E (MO  . 1 i;.)OE  - i 90. 1 OOOf  - 1 9 0 . 8 « 7 5E-U  1 U . 2 ->4  J t 0 0 0 . 1 0 0 0 c.  - 1 9 
0 . .3  2 3 2k.  0 0 0 . 2 4 i nr,  n 4 0 . 1 0 0 3 E - 1 9 o . ] uOOE-  1 90 . 9 1 5 5E - U 1 0 . 2 3a  Sr.  000 . 10UOE- 1 9 
0 .( 1 6 2 k.  000.2  1 7-..-'  0“') . I OO.iE  - j 70. 1 OOOb-1  90. 88  3 1 k.-Ol  0. 1 6641  000 . 1 OOOE  - 1 9 
0 . 54  1 9b  0 0 0.1  452 r.  0 *. 0 . 1 0 0 0 E - I 9 0 . 1 0 0 0 F - 1 9 0 . 3 2 b 3 E - 0 1 0 . U'  1 5 1.  000 . 1 OOOE- l 9 
0 . / 12  3 u.  0 o 0 . *>  1 o < .*  O 3o  , I 0:).)“,  - J 40 . 1 book  - t oo  . 1 o7  2r.-0  10  . j9  / 6c  -0  1 0 , 1 uOOk,  - 1 4 
0.901  3K  OoO  . 5 i n7r.  .)  jo  . i 0()9r  - l 9 o . 1 00  Ok.  - 1 90 . 7 o5(ii:.-o20  . J 639b  - 0 l U . 1 00  Ob  - 1 9 
o . I )"(!'■.  >10.  O , 1 'H'i.i  t - i 4 o . 1 O II o i-  - I V 0 , J 3 9 - o 2 0 , 7 M> f -020, 1 o 0 0 b - 1 0 

o . 1 2 5 4 f,  I)  1 o . 3 j / 3).  o 7 1) . 1 00  )E  - 1 70 , 1 00O1  - t ^ J . I50.lr.-0  2o.  35  3 95-070  . loOOK-  1 9 
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0.1423E  0 1 0 . 3 1 ObE 
0.1579E  0 1 0 . 2 9 b 7 E 
0.1732E  010. 2 904 E 
0.27  UK  01 
0.0  0.2527E 

0. 8041 E -010. 2 5 2 6 E 
0.160HE  000. 252 2 E 
0.2415E  0 0 0 . '2  5 0 7 K 
0.3250E  000. 244 3 E 
0.4201E  000.2176E 
0.5459E  000 . 1 53 1 E 
0.7109E  000. 919 9 E 
0.8964E  000.5924E 
0. 108215  010.44 0 3 E 
0.1258E  0 1 0 . 3 o b 4 E 
0.1  425E  010. 329 0E 
0.1584E  010. 3086 E 
0.1738E  010. 2 » 72 E 
0.1888E  010,2912 E 
0.318E  01 
0.0  0.2533E 

0.8049E-010.2532E 
0.1610E  0 0 0/ . 2 5 2 7 E 
0.2422E  000.2509E 
0.3271E  000.2440E 
0.4244E  000. 2 178 E 
0.5498b  000.1597E 
0.7109E  0U0.1013E 
0.8933E  0 0 0 . 6 b 4 3 E 
0.1078E  010. 4685 E 
0.1256E  010. ■ 3982 E 
0.1423E  0 10. 3 5 0 4 E 
0.1587E  010. 323 2 E 
0.1  743E  010. 30 67 E 
0.1895E  010.2967E 
0.2043E  010.291  IK 
0.358E  01 
0.0  0.2539E 

0. 8056E-01 1) . 253  7 E 
0 . I b 1 2 E 0U0.2532E 
0 . 2 4 3 1 E 0 0 0 . 2 5 1 2 E 
0.3295K  000. 2437 K 
0.4288E  000.21  ?8E 
0 . 5 6 3 7 E 000  . 1649E 
0.7  U3E  000.109 5 E 
0.8905E  0O0.7337E 
O.  1 075E  010.53.-j1E 
0.1264E  01O. 4334b 
0.  U26E  010.374  IE 
0.1  590E  0 1 0 , 3 3 5 9 E 
0.1748E  010, 31 8 oE 
0 . 1 9 0 1 E 0 1 0 . 3 0 4 j E 
O.2091  E 0 1 0 . 2 9 o 1 E 
0.2198 E 01 0,2 90 9K 
0,398c.  0 1 
0.0  0 . 2 5 4 t E 

0.8064E-010. 2543k. 
0 , 1 6 1 o E 0 0 0 , 2 5 3 7 1-. 
0 . 2 4 4 2 E 0 0 0 . 2 5 1 5 E 
0.3  32  2E  0 0 0 . 2 4 3 J E 
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0.2384c  010.3285c  030. I OOOE-1 90. 10 00 E-19 D.1319E -020. 218 4 E-0 20. 1000 E-19 
0.2532c  01 0 . 31 22 E 0 30 . 1 000c- 1 90 . 1 0 uOt- 1 90 . 8 9 50 E-0 30 . i <W4b-020 . 1 OOOE- 1 9 
0 . 26  7 9E  010.303  6b  0 9 0 , l OOOE- 1 4 0 . 1 00 OH  - 1 90 . 5 7 3 U.-0  90 . 9b2 2E-0  3 0 . 1 OOOE- 1 9 
0.28  2 3c  0 10. 20 7 OK  0 3 0 . loOOE- 1 9 0 . 1 OOOE- 1 90 . 3 3 1 2k -0  30 . 558 3 E-0 30 . 1 0 00c-  1 9 
0 . 2 9 8 o c 01  0.2920b  0 30 . 1 OOOE-1  9 0 . 1 OOOE-  1 9 0 . 1 48  .iE-0  30 . 2 805E-030 . 1 0 0 OE-  1 9 
0.644E01 

0.0  0.2  58  IE  040 . 1 OOOE-l 90. 1 OOOE-1 90. 9*4BE-0 I 0. 2424E  00 0 . 1 00 OE- 1 9 

0 , d 3 3 9 E - 0 1 0 . 2 5 7 7 1 04 0 . 1 OOOE- 1 90 . 1 OOOK - 1 90 . 1 002c  00O.2J98E  000.1000E-19 
0. 168  IK  OOO.  2558k  040. 1000E-1 90  . 1 00ftL-19O. 1 J22E  00O.2313E  OOO.lUOOK-19 
0.25/4E  000. 2 8 03 E 040 . 1 OOOE- 1 90 . 1 OoUb- t 90 . 1 04 3c  00O.2150E  O00.1000E-19 
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0.3539c  000. 2375K 
0.402215  000. 2 J4 -lfc, 
0.58566  0 0 0, 18  16 E 
0.72745  000.14526 
0 . 8 H 6 6 6 0 0 0. 1 1 2 56 
0.10586  01  0,8747b 
0.1  234c  010. 6 9935 
0.14105  010.58015 
0 . 1 5 b 4 c 0 10.4988b 
0.1753E  010.44215 
0. 1918E  010. 4015 K 
0.20786  010,37186 
0.2235b  010.3499b 
0.23875  010.3332b 
0.2537b  010.3202b 
0.2685b  010.3099b 
0.28305  010.30175 
0.2973E  010.2955b 
0.31165  010.29115 
0.6856  01 
0.0  0,2584b 

0.8430E-010.2579E 
0.1700E  000.25555 
0.2604E  000.2493E 
0.35815  000.23605 
0 , 4 ft 6 9 E 000.21345 
0.5907E  000.18235 
0.7314E  000.14795 
0.8882E  000.1164b 
0.1057c  010. 9151 E 
0.1232E  010.73545 
0.14085  010.6104E 
0.1581E  010.5238b 
0. 1752E  010,46275 
0.1917E  010.4185b 
0.20795  010.38585 
0,2236c  010.3ol4E 
0.2390E  010.34285 
0.25425  010.32835 
0.2690E  010.31695 
0.28365  U1 0.30775 
0.2981E  010.3002b 
0.31245  010. 2 9 37 E 
0 . b V 4 5 01 
0.0  0.2585G 

0.86466-010. 25695 
0.1731c  000.25645 
0 . 2 6 4 8 £ 000. 2 3*0 r, 
0.35945  000.2357b 
0.59265  000,1824b 
0.8878E  0 0 0 . 1 1 7 2 K 
0.1232b  010.74805 
0. 1681E  010.52925 
0.191/5  010.42225 
0.22375  010. 3ft 3 9b 
0.25435  010.3301b 
0.243/5  u 1 0 . 30  6 8 k 
0.3125c  0 1 0 , 2 9 4 1 5 
0.7415  (H 
0 . 0 0.25811; 
0 . 8 o 4 ft  5 - 0 1 0 . 2 5 o 0 5 


040. 100 05-  l 90. 1000b- 1 90. 10 39E  000. 189 4E  000.10005-19 
040. 1000E-I90.10u0b-I90. 97155-010. 1557b  000.1000E-19 
040. 10 07E-190. 100 06 -19 0.B119E-010. 11896  000.10005-19 
04  0. 1O00E- 190. 1000  b- 19 0.59 55 5-0 10. 8 4555-0 10. 1000E-  19 
0 4 0. 10  0 OS-190. 1000 5-190. 3996 K -010. 5 6 80t-010. 10005-19 
030. 1 0 0 0 b - 1 9 0 . 10  005-19  0.2  616  5 -010, 375  25 -010. 10005-1  9 
0 30. 10005-1 90. 1000b- 190.  17405-010. 25315-010. 10005-19 

030. 10005- 1  90  . 10  005-190. 1 1905-010. 1761 t-Q 10. 100  05-19 

030. 100 0 5-1  90. 100 05-190. 83b7E-020, 12585 -010. 10055-19 
0 30. I 0005- 190. 10005- 190. 60065-020. 9 1765-020. 100  Ob- 19 

030. 10005- 1 90 .1000b- 190 .43745-020.67785-020 . 10005-19 

030. 10005- 1 90.10005-190. 32085-020. 50365-020. 10005-19 

03 0.100 05- 190. 10005- 190. 23525-020. 3733E-020.1000E-19 
0 30 . 10  005- 1 9 0 . 10005-190 . 17  0 8 5-020. 2 7 37 5-0  20. 10005-19 

030.10076- 190.10005-190.12125-020.19575-020.10005-19 
030 .100  05-1 90. 1 0 005-1  VO. b 2 095-030. 13345-020. 1 0005- 19 
030. 1000-6-1  90. 10  005-  190. 51  17  5-0  30. 8 3 525-030. 1000b-  19 
0 30. 10006-190. 10 OOE. -190. 2782 5 -030. 45545-030. 1000 5 -19 

0 30. 10006- 1 90. 10005- 19 0.1 157 6- 030. 169 76-030. 1000 E- 19 

040. 10005- 190. 10006-190. 1015E  000.23786  000.10006-19 

040.10006- 190.10006-190.10226  000.23505  000.10006-19 

040.10006- 190.10006-190.10396  000.2260E  000.1000E-19 

040.10005- 190.10006-190.10545  000.2095E  000.1000E-19 

040. 10006- 1 90. 10005- 190. 1045E  000.18456-000.10005-19 

0 40. 1 0006-  190. 10 006- 190. 97 395-0 10.1 5 25E  000.10006-19 
0 40. 100 06-  1 90. 10  005-190. 8 2 076-01 0. 1 1 806  000.10006-19 

040. 1000 6 - 190. 100 06-190. 6151 6-0 10. 8 5856-01 0.10006 -19 

040.10006- 190.10006-190.42395-010.59246-010.10006-19 
030. t 0006- 1 90. 10006- 1 90.28385-010.40016-010. 10006-19 
030. lOOOE-l 90. 100 0h-190. 19 206-010. 2741 6 -010. 10006-19 

030. 10005- 190.10006-190.13326-010.19305-010. 10006-19 

030. 10006- 1  90.  10006-1  90. 9 4 856-020. 1394E-010. 1000 6 -19 
03 0.40  0 06-190. IOuO 6-190. 68 9 3 6-020. 10286-010. 1000 1-1  9 

030 .10006-  190. 1000b- 190. 508 55-0 20. 7 68 16-02 0.1 0006- 19 

030. 10 00 6- 1 90. 10006-190. 3 7 83 6-020.5784 6-020. 1000 6-19 
0 30. 100.76-1  90. 10006-190.2  3 226-020. 4 3606-020 . 100 OK-  19 
03 0. 3 0 0 Ob-1  90. 10  00  6- 190. 20966-020. 326  85-0  20. 10  0 0 E- 1 9 

0 3 0. 1 0076- 1  90. 1000 6- 1 9 0.15376-020. 24156 -02 0. 1000 E -19 

030 . 10006-  1 90. 1 0006-190. 10996-020. 1 7376-020 . 10006-1 9 
0 30. 1 00  06- 1 90. 10006-1 90.7 4826-030. 11 895-020. 10005-19 

030. 10006- 190. 10006- 190. 46066-030. 73456-030. 10005-19 

030. 10006- 190.10006-190.21666-030.34606-030.10006-19 

040.10076- 190.10006-190.10206  000.23676  000.10006-19 

040 . 10076- 1 90 . 1000K-1 90 . 1 0316  000.23006  000.10005-19 
040 . 1 90 76-1 90 . 10006- 1 90. 10425  000.2249E  000.10076-19 
04u.. 1 0005-1  90. 10006-190. 10446  000.20006  000.1  0006-19 

040.10005- 190.10006-190.10466  000.18346  000.10006-19 
04  0. 10  00b- 1 90. 10006- 190. 822ft 6-01 0. 1 17dt  00  0.10006-19 
040. 1000E-1 90. 10006-190. 42906-010.59726-010. 1000E- 19 
0 10. 1 00  0 b-  1 40. 10006-190. 19  596-61  0.27  8 bf.-O  10. 10006-1  9 
030. 1 00  06-1 90. 10  906-190. 9 7 28  6-0  20. 14236-01  0. 3 00 OE-1 9 
0 39. 1000  6- 1 90. 1 0U06- 190. 62405-020, 7 8 7 7 6-020. 1 0006-  19 
0 3 0 . 1 0 0,7  6-  I 90. 1 000b- 190. 2 927  6-020. 4 4986-0  20. 100  06-19 

030. 10006-  190. 10006- 190.160/6-020. 25 125-020. 10006- 19 

030. 1 0076- 1 90. 10008 -1 90.78936-030. 12476-020. 10005-19 
0 3 0 . 1 0 0 0 6 - 1 9 0 . 1 0 0 0 H - 1 9 0 . 2 2 9 7 6 - 0 3 u . 3 6 4 8 c - 0 3 0 . 1 0 0 0 £ - 1 9 


040. ] 0076-1 90. 10006- 190. 1 0406  000.23096  000.10006-19 
040. 100.76-1  90. 1 0006- 190. 10486  0 0 0. 22  306  000.100  06-19 
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0.17546  000.25436  040 . 1 OOOE-1  90  . t OOOE-  1 90 . 1 0 7 hi  000.21866  000. 10006-19 

0.25966  000.244(31-:  040 . 1 OOOE-  1 90  . 1 0u0fc>l  9 0 . 1 0506  000.19806  000 . 1 00  0 6>  1 9 

0.36396  000.23396  040 . 1 00  06-  1 90 . 1 0006-  1 90 , 1 04  96  000.17801-3  000.10006-19 

0.59616  000.1027 6 0 40. IO00E-190. 10  0 0b-190, 8317 E -010. 1161 R.  000.1  000E-19 

0.89096  000.  120  bE  0<t0 . 1 0006-  1 90 . 1 0006-1 90 . 44946-0  1 0 . 6 1 / 36-0 1 0 . 1 00  OK-  1 9 
0.12296  010.78126  030 . 1 00 OE- 1 90 . 1 000 6- 1 90 . 2 1 6 1 6-0 1 0 . 30 1 66-0 l 0 . 1 0006- 1 9 
0.15766  010.55806  0 30 . 1 0006- 1 9 0 . 1 0006- 1 90 . 1 1 086-0 1 0 . 1 5 866-0 1 0 . 1 00 06- J 9 
0. 191 4E  010.44206  030.10006-190.10006-190,61026-020.89516-020.10006-19 
0.22366  010.37756  030. 10006-190. 10006- 190. 34926- 020. 52316-0 20.10006-19 
0.25456  010 . 33 9 8 6 030.10006-190. 1000 6 -190. 2 0026-020. 30486-02 0.10006-19 
0.28436  010.31686  030. 10006-190.10006-190. 11076- 020. 170 56-02 0.10006-19 
0.313  36  010.30276  0 30 . I 00 06- 1 90 . 1 0006- 1 90 . 56 1 76-0 30 . 8 7 3 36-0 30 . 1 00 06- 1 9 
0.34196  010,29406  030.10006-190.10006-190.22736-030.35536-030.10006-19 
0.7856  01 

0.0  0.25766  040. 1000E-190. 10006-190. 10596  000.22516  000.1000E-19 

0. 88826-010. 25546  0 4 0 . J 00 06- 1 90 . 1 OOOfc - 1 90 . 1 0656  000. 21806  000.10006-19 
0.17766  000.2532?:  0 40 . 1 0 0 0 6 - 1 90 . 1 0006-  1 9b  . 1 0 7 1 6 000.21266  000.10006-1  9 
0.26726  000.24306  0 40 . 1 00 06- 1 90 . 1 0006-1 90 . 1 0606  000.19206  000.10006-19 
0. 36826-000. 2324E  04 0 . 1 00 06- 1 90 . 1 0006- 1 90 . 1 0 53 6 000.17306  000.10006-19 
0.60356  000.18266  040.10006-190.10006-190.83766-010.11456  000.10006-19 
0.89416  000.12326  04 0.10006-190. 10006-1 90. 46846-010.63256-010.10006-19 
0.122  96  0 10. 8 1356  0 30 . 1 0 0 06- l 9 0 . 1 0 006- 1 9 0 . 2 3 3 46-0 1 0 . 320 76-0 1 0 . 1 0006-1 9 
0.15746  010.58336  030.10006-190.10006-190.12286-010.17256-010,10006-19 
0.19126  010.46026  030. 1 00 06-1 90. 10006- 190. o 9056-020. 9932 6-0 20. 10006-19 
0.22366  010.39066  030.10006-190.10006-190.40526-020.59476-020.10006-19 
0.2547  6 010.34976  0 3 0.10006-1 9 0. J 0006-190. 24 136-0  2 0. 35996-02 0.1000E-19 
0.2847  6 010.324  46  030 . 1 0006- 1 90 . 1000E-190 . 1 4 1 2E-020 . 2 1 326-020 . 1 0006- 1 9 
0.31406  010.30826  030. 1 0006-1 90. lGOot-190. 7800 6-030 . 11896- 02 0.10006-19 
0.342  76  0 10.29606  0 3 0 . I 00 06- 1 90 , 1 0006- 1 90 . 3 8 57 6-0 30 . 592 1 6-0 30 . 1 0 006-  1 9 
0.37106  010.29196  030.10006-190.10006-190.14996-030.23106-030.10006-19 
0.8306  01 

0.0  0.25586  040.1 0006-1  90. 1 0006-3  90. 10776  000.21906  000.10006-19 

0.89806-010.25446  040 . 1 0006- 1 90 . 1 0006-1 90 . 1 0806  000.21206  000.10006-19 
0.17986  000.25206  040 . 1 00 06- 1 90 . 1 0006- 1 90 . 1 0856  000.20656  000.10006-19 
0.27006  000.24146  040 . 1 0006- 1 90 . 1 0006- 1 y 0 . 1 0 / 06  000.18706  000.10006-19 
0.37266  000.23086  040 . 1 00 36- 1 90 . 1 0006- 1 9 0 . 1 0 566  000.16816  000.10036-19 
0.60916  000.18286  040 . 1 0006- 1 90 . 1 0006- 1 90 . e 4 1 66-0 1 0 . 1 1 296  000.10006-19 
0.89756  000.12576  040.10006-190.10006-190.48586-010.64406-010.10006-19 
0.12  296  010.84296  030 , 1 0006- 1 90 . 1 00  06- 1 90 . 2 4 966-0 1 0 . 3 3 7 46-0 1 0 . 1 0 0 06-  1 9 
0.15726  01  0 . 6075  K 0 3 0 . 1 00  :)S  - 1 90 . 1 0006-  1 90  . I 34  26-0  1 0 . 1 8 556-01  0 . 1 00  36-  1 9 
0.1  91  06  010.477  76  0 30,  t 5006-190. 10006-190. 7 6916-020. 108  7 e.-OlO. 10006-19 
0.22356  010.40396  03 0.10 006-190. 10006-19 0.4 6106-02 0.56 4 66-0 2 0.1000 6-19 
0.25496  010.35946  030 . 1 0006- 1 90 . 1 000 t - 1 90 . 28 206-02 0 . 4 1 326-020 . 1 0006-  1 9 
0.28516  010.33196  030 . 1 0006- 1 90 . 1 0006 - 1 90 . 1 7206-020 . 2 5526-020 . 1 0006- 1 9 
0. J J 4o6  010.31416  030. 10006-]  90 . 10006-1  90. 10176-020. 15256-020. 10006-19 
0. 34356  010.30246  030 . 1 0006- 1 90 . 1 0006- 1 90 . 55696-0 30 . 84 l 1 6 -0 30 . 1 0006- 1 9 
0.3  7196  010.29426  030 . 1 0006- 1 90 . 1 0005- 1 90 . 24096-030 . 36546-0 30 . 1 00 06-  1 9 
0.8766  01 

0.0  0.25596  040. 10006-190. 10006-190. 10936  000.21286  000.10006-19 

0,91  186-010. 25306  040 . 1 0006’-  1 90  , H>006-  1 90 . 1 0 956  000.20706  000.10006-19' 
0.18216  000.25086  040. 10006-1 90. 1 0006-190. 10976  000.20056  000.10006-19 
0,2777c.  00(1.24006  040 . 1 00  36- I 90 . 1 0006-  1 40  . l 07  06  000.18206  000.10036-19 
0.37706  000.2291  r..  040.1  0006-1  40.10(106-190.10586  OOO.loJSE  OOU.  10006-19 
0.6143c,  000.1  8276  040. 100  )6-1  90. 1 0 0 0 5 - 1 9 0 . 8 4 q 5 6 - 0 1 0 . 1 1 1 36  000.10006-19 
0.901  16  00  0.1  277 6 040 . 1 000 1- 1 90 . 1 000 1 - J 90 . 50 1 16-0 1 0 . 65256-0 1 0 . 1U0 06- 1 9 
0. 1 2296  010. 86926  030. 10006-1 90. 1 0006- 1 90 . 26476-5 1 0 . 35 1 96-0 1 0 . 1 0006- 1 9 
0,15  70c.  010.62925  030, 1 6 506-19 O.lOuOl-19 0.1 45 16- 015. 19 756-01 0.10006-19 
0,1909c  0 1 0 . 49486  0 39 . 1 oO 06- 1 90 . 1 0u06- l 9 0 . 5 4 9 7 6-0 2 0 . 11 7 / 6 -0 1 0 . 1 0 0 OE- 1 9 
0.22386  0 10.41896  085. 1 HOOK  - 1 4O  . 1 0006-1  90 . 5 1 91)6-020 . 73166-020 . 1 0006  - 1 9 
0.255  1 6 010  . 4 o 9 2 6 0 30.1  0036-  1 90  . 10001  -1  40  . 3 2226-0  20  . 4 64  56-0  20  . 10036-19 
0 . 2 8 5o  6 0 10.33921,  030.  I 0006  - l 90. 1050c.-  1 90 . 20  l or.-O  20 . 254  86-020 . 10006-19 


0.31526  01  0.3  mu 
0.34416  010.30616 
0.37276  010,29726 
0.40096  010.2917E 
0.92  IE  01 
0.0  0.25  4 8 E 

0.91805-010.25206 
0.18445  000,24936 
0.28245  000.23806 
0.38145  000.22755 
0.61945  000.1825E 
0.90515  000.1295E 
0.1230E  010. 8932 E 
0.1570E  010.650UE 
0.1908E  010.51105 
0.22365  010.42925 
0.25535  0 1 0 . 3 7 8 6 E 
0.2859E  010.3463E 
0.31576  010.3250E 
0.3448E  010.31065 
0.37355  010.3007E 
0.4016E  010.29386 
0.3665  01 
0.0  0.25355 

0. 92945-010. 25005 
0.18685  000.24775 
0.28615  000.23665 
0.33565  000.22575 
0.62496  0 00. 1 3 2 1 E 
0.90905  000. 1310 E 
0.12325  0L0.9150E 
0.1559E  010.6696c, 
0 , i 9 0 7 E 010.52645 
0.2235E  010.4410c 
0.2554E  010.38775 
0.2861E  010.35345 
O.ilblE  010.33055 
0.34535  010, 3 14 7 E 
0.3741  E 010.303 4 6 
0.40245  010.29495 
0.4303E  010.29016 
0.1015  02 
0.0  0.25216 

0.95UE-010. 24906 
0 . 1 8 9 1 S 0 0 0.24605 
0.26615  000.2350  6 
0.39015  000.22396 
0.63005  0 0 0. 1 6 1 1> 5 
0.91295  000.13235 
0.12335  010.93456 
0.15655  010.607hr 
0 . 1 9 0 o c o 1 0 , 54  096 
0.22355  010.45235 
0.28645  0 10.390  4 c. 
0.2  1635  0 l o.3  6i'2k 
0 . 3 1 o 4 6 0 1 0 . 3 3 1 7 6 
0.345Hc  0 1 0 , 3 i b 9 5 
0 , 3 7 4 7 c u 1 0 . 3 0 h 9 k 
0.4  9325  01  0.298  o I' 
0,4311c.  0 l i) . 2 9 2 ii  c 
0.1085  02 


030. 10005-1  90  . 1 0006-19 0.1 2 32 5-020. 1 81  95-020, 1 000E-19 

030. 10005- 1  90. IOuOE- 190. 7 0695-030. 1051 E-020. 1000E- 19 
030. 1000E-1  90. 10005-190.36055 -030. 5390 5 -030. 10005-19 
030. 1000E-190. 1000 5-190. 1433 5-030, 21505-0 3 0.1000E-19 

040. 1000E-1 90. 1000E-190. 1 108E  000.20665  000.1000E-19 

040.10005- 190.10005-190.11075  000. 20055  000.1000E-19 
040. 1000E-1 90. 1000E-190. 1 107E  000.19455  000.1000E-19 

040 . 10005- 1 90. 10006-190. 10835  000,17675  000.10005-19 
040. 1000E-190.100U6-190. 10596  000.1589E  000.10006-19 
040. 1000 E -190. 10006-190. 8 460 E- 010. 10966  000.10006-19 
040.1 000E-1 90. 10006-1 90. 5 1436-0 10. 656 46-01 0. 10005-19 
030. 1000E-1 90. 10006-190.27876-01 0.3 6486-010, 10005-19 

030.10005- 190. 10006-190. 15 56E- 010. 2085E- 010. 1000 E- 19 

030. 1 0005-  190. 10006- 190. 9 201 c-020. 12615-0 10. 10006-19 

03 0.10006- 1 90. 100 06- 19 0.56 9 45-020. 7 955 6- 020.10005-19 
030. 1003E-190. 1000 6-190. 36155-02 0.51355-020. 10005-19 

030.10005- 190.10006-190.23125-020.33286-020.10006-19 
030. 100 OE-1 9 0.100 06-1 9 0.14616-020. 21 2 55-020. 1000E- 19 

030.10006- 190.10006-190.88966-030.13056-020.10006-19 

030. 10005- 190,10006- 190. 49815-030.73526-030.10005-19 

030.10005- 190. 1000 E- 190. 218y 5-030. 3242 5-030. 1000E-19 

040.10006- 190.10006-190.11205  000.20036  000.10005-19 

040.10005- 190.10005-190.11185  000.19456  000.1000E-19 

040. 10006- 190. 10006-190. 11166  000.18865  000.10005-19 

040.10005- 190.10005-190.10876  000.17106  000.10006-19 

040. 10005- 190. 10006-190. 1089E  000.1545E  000.10006-19 

040.10006- 190.10006-190.84625-010.10796  000.10006-19 

040.10006- 190.10006-190.52576-010.66236-010.10006-19 

030.10005- 190.10006-190.29176-010.37606-010.10006-19 

030. 1000 6- 190. 1000 6-190. 1656 6-0 10. 2 lb6E-010. 1000 E- 19 

030.10006- 190.10005-190.99246-020.13406-010.10006-19 
030. 1000E-1 90. 100 06- 190. 62 2 05-020. 85666-0 2 0.10006-19 
0 30.  10  006-1  90. 10006-190.40066-020.561  16-020. 10006-19 
03 0. 10 0 0 £-190. lOOOt-l  9 0.2 609 E-02  0. 37046-02 0.10006-1  9 
030 . 1 0006  - 1 90'.  10006-19  0. 16876-0  20. 242  26-020 . 10006-19 

030. 10006- 1 90; 10006-190. 10566-020. 15296-020. 10006-19 

030. 100 06 - 199. 10006-190. 608 26 -030. 88616-0 30. 10006-19 

030. 10006-  1 90. 10006-190. 27 30c- 030. 399 16-030 . 10006-19 

030.10006- 190.10006-190.84376-040.12366-030.10006-19 

040. 1 0006- 190. 10006- ) 90.1 1316  000.19426  000.10006-19 
040. 10036-190. 1000E-190. 11276  000.18906  000.10006-19 

040. 10006- 190. 10006-190. 1122E  000.18286  000.10006-19 
0 4 0 . 1 0 0 0 6 - 1 -9  0 . 1 0 0 0 1,  - 1 9 o , 1 U y 0 6 00  0.1  66  06  000.10006-19 

040.10006- 190.10006-190.10576  000.15026  000.10006-19 
040. 19006-1 90. 10006-190.84506-010.10616  000.10006-19 

040.10006- 190. 10 006- 19 0.53 5 16- 01 0.6 64 3c-0 10. 10006-19 

030. 10006- 190. 1000  £-1  9 0,3  0 355-010. 38576-01  0.1000  6' -19 
0 30. 1 00  0E-  1 90. 1 000H -1  90. 17 506-0 10. 227 86-0 10. 1 0006-1  9 
030. 1 OOOc-l 90. 10006-190 . 1 062 K-0 10. 1 4146-010. 10006-1 9 
0 3 0 . 1 0 0 0 6 - 1 9 0 . 1 0 0 0 6 - 1 9 0 . 6 7 3 0 6 - 0 2 0 . 9 1 4 4 6 - 0 2 0 . I 0 0 0 6 - 1 9 
0 30 . 1 00 Oc- 1 90. 10006-1  90. 438  76-020.606  36-020 . 10006  - 19 

030 . 10006-  1 90. 10006-1 90. 28976-020. 40606-020 . 10006-1 9 
03 U . 10006- 1 90. 1 0006-190. 1 9 0 4 6- 0 2 0 . 26 996 - 0 2 0 . 1 0006-19 
0 8 0 . 1 0 0 i)  c.  - 1 9 0 . 1 u 0 0 6 - 1 9 0 , 1 2 2 3 6 - 0 2 0 . 1 7 4 9 6 - 0 2 0 . 1 0 0 0 6 - 1 9 
0 3 0 . l on  Or  - 1 90  . 1 000  c - 190 . 7 4 886-030 . 1 07«6-O20  . 1 0006  - 1 9 
0 3 0 . 1 00  06-1  90  . 1 0006-  190. 4 1 7 o 6 - 0 3 0 . 6 04 2 6 -0 3 0 . 100  06-19 
0 30 . 1 0 0 0 6-]  90. 1 0i)0 c.-l  90. 1 82  36-030. 26  456-030. 10006-  19 
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3 


i 


0.0  0.25U46  040.  in00E-J9U.10O0E-190.H39e.  000. 1 8806  000.10006-19 

0.94326-010.24706  040 . 1 0006-  1 90 . 1 OOOfc  - 1 90  . II  3 26  000.18256  000.10006-19 
0.19146  000.24416  040 . 1 00 96- 1 90 . 1 0006- 1 90 . 11 266  000.17706  000.10006-19 
0.29286  000.28306  040 . 1 00 06- 1 90 . 1 0006- 1 90 . 1 0906  000.16156  000.1000E-19 
0.39466  000.22206  040 . 1 00060-  1 90 . 1 0006-  1 90 . 1 06360  000.14606  000.10006-19 
0.63986  000.16016  0 40 . 1 0006- 1 90 . 1 0006- 1 90 . 8 j9 1 6-0 1 0 . 1 0 256  000.10006-19 
0,92086  000.1341  6 040.10006-190.  1 OuOE- 1 90 . 54 9 2 6-0 1 0 . bn 396- 0 1 0 . 1 0006- 1 9 
0.12376  010,96806  0 30 . 1 00 06- 1 90 . 1 0006- 1 90 . 3 24 26-01 0 . 40 1 26-0 1 0 . 1 00 OE- 1 9 
0.16706  010.72066  030.10006-190.10006-190.19266-010.24406-010.10006-19 
0.19056  010.56826  03 0.10 U0E-1 9 0.1000 6-100.11956-010. 15 506-010. 10006-19 
0.2  2 356  010.47416.  030 . 1 00 06-  1 90 . 1 0006.-  1 90 . 7 7 246-020 . 1 02 36.-0  1 0 . 1 0 006-  1 9 
0.2556  6 010. 413 9 E 030. 10  07  6-190. 1 0 00 6-19 0.51  4 06 -020. 69 2 46-020. 1000 6 -19 
0.28676  010.37376  030.10006-190.10006-190.34736-020.47496-020.10006-19 
0.31706  010.34646  030. 100 06-1 90. 10006-1 90. 23546-020. 3 25 76-020. 1 0006-1 9 
0.34676  010.32726  030. 10006-1 90.10006- 190. 157 46-020. 2 2006-0 20.10006-19 
0.3757E  010.31346  030. 10 OOE -190. 100 06- 190. 10166-020. 14286-020. 10006-19 
0.40446  010.30306  030.10006-190. 100 06-190. 598 06- 030. 84666-030. 10006-19 
0.4327  6 010.29486  030. 100  06.-190. 10006- 190. 2 72  7 6-030. 387  26-030.10006- 19 
0.46066  010.29016  030. 1000 6-190. 10U06- 190. 8526 6 -040. 12 12 6- 030. 10006- 19 
0.1106  02 

0.0  0.2435E  040. 10006-190. 1 0006-1 90 . 11 4SE  000.18206  000.10006-19 

0.96286-010.24536  040.10006-190.10006-190.11366  000.17676  000.10006-19 
0.19386  000.24216  040.10006-190.10006-190.11286  000.1714E  000.1000E-19 
0.29596  000.23106  040.10006-190.10006-190.10906  000.15706  000.10046-19 
0.39896  000.22006  040.10006-190.10006-190.10496  000.14196  000.10006-19 
0.63516  000,18096  040.10006-190.10006-190.84266-010.10436  000.1000E-19 
0.91696  000.13336  040.10006-190.10006-190.54306-010.66486-010.10006-19 
0.12356  010.95226  030. 10 00 6 -190. 10 006- 190. 31436-010. 3 9416-0 10. 10006-19 
0. 15696  010 ,704bt£  030 . 1 0006-  1 90 . 1 0006- 1 90 . 1 84  06-0 1 0 . 2 3636-0  1 0 . 1 0006- 1 9 
0.19066  010.58496  030. 1090 6-1 90. 1 0 006- 1 90. 1 1296-01 0,1 48 46-010. 1000E- 19 
0.2235  6 010.4634  6 0 3 0. 10006-190. 10  006- 1 9 0. 7232  E- 02  0.96  9 7t,-020. 100  06-  19 
0.2  8 556  010.40536  0 30 . 10 0 96-1  90. 10006- 190. 4 7 66 6- 020. 65026-02 0.1 0006-1  9 
0.28656  010.36696  030. 100  06-190. J 0 006-1 90. 3186 6- 020. 44096-020. 10006-19 
0.31686  010.3411E  030 . 1 00 OE - 1 90 . 1 0006- 1 90 . 2 1 306-0 20 . 29836-020 . 1 00 OE- l 9 
0 . 34626  010.32316  030 . 1 0006-1 90 . 1 OuOc  - 1 90 . I 4026-020 . 1 9b 26-020 . 1 0006- 1 9 
0.37526  010.31036  030. 100 06- 190. 10006-1 9 0.08816-03 0.128 46 -020. 10006-19 
0.40  366  010.30  106  030 . i 0 006-1 90 . 1 0006- 1 90 . 5 1 52h~03o  . 7 37 36-030 . 1 OOOE- 1 9 
0.4319  6 010.2  9396  0 30 . 1 0 096-1  90 . 1 Ooot-  1 9 0 . 23236-0  J()  . 33356-0  30 . 1 0006-1 9 
0.0 
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NWAVl  2000  2500 


►V  A V 

2000.0 

0.0 

0.0 

0.0 

0.13b- 

03 

0.0 

0.13E 

00 

0.66E 

00 

WAV 

2000.0 

0.25E 

01 

0. 1 5E 

0 2 

0.185 

03 

0 . 4 0 E 

0 3 

0.37E 

03 

0 . 2 8 E 

04 

0 . 7 0 E 

04 

WAV 

2025.0 

0.0 

0.0 

0.0 

0.78E- 

03 

0.0 

0.26E 

00 

0. 1 IE 

01 

WAV 

2025.0 

0.25E 

01 

0.  15E 

02 

0 . 1 8 E 

03 

0.40E 

03 

0 . « 7 E 

03 

0 . 2 8 K 

04 

0.63E 

04 

WAV 

2050.0 

0.0 

0.0 

0.0 

0.37b- 

02 

0.0 

0. 50 E 

00 

0 . 1 7 E 

01 

WAV 

2050.0 

0. 2 5E 

01 

0.  15E 

0 2 

0 . 1 8 5 

03 

0.40E 

0 3 

0 . 8 7 E 

03 

0.27  F. 

0 4 

0.51E 

04 

WAV 

2075.0 

0.0 

0.0 

0.0 

0.15b- 

01 

0.0 

0.95E 

0 0 

0.26E 

01 

WAV 

2075.0 

0.25E 

01 

0 . 15E 

02 

O.iriE 

03 

0 . 4 0 E 

0 3 

0.87E 

03 

0 . 2 7 E 

04 

0.44E 

04 

WAV 

2100.0 

0.0 

0.0 

0.0 

0.52E- 

01 

0.29E 

00 

0 . 1 7 E 

01 

0.37E 

01 

WAV 

. 2100,0 

0.25E 

01 

0.1  5E 

02 

0.185 

CT3 

0.4  0E 

0 3 

0 . 86E 

03 

0.25b 

04 

0.36E 

04 

WAV 

2125.0 

0 . 0 

0.0 

0.0 

0.18b 

0 0 

0.69E 

00 

0.28E 

01 

0.50E 

01 

WAV 

2125.0 

0.25E 

01 

0.  15E 

02 

O.tdS 

03 

0 . 40E 

03 

0 • 8 2 E 

03 

0.21E 

04 

0.29E 

04 

WAV 

2150.0 

0.0 

0.0 

0.0 

0 . 5 2 E 

00 

0. 15E 

01 

0.45E 

01 

0 . 7 3 E 

01 

WAV 

2150.0 

0 „ 2 5 E 

01 

0 . 1 5 E 

02 

0.185 

03 

o 

• 

o 

•TJ 

03 

0.75E 

03 

0 . 1 7 E 

04 

0.23E 

04 

WAV 

2175.0 

0.0 

0.0 

0 . 3 1 E 

00 

0 . 14E 

01 

0 . 3 1 E 

01 

0.67E 

01 

0.87E 

01 

WAV 

2175.0 

0.25E 

01 

0. 15E 

0 2 

0.185 

03 

0.395 

03 

0 . 5 5 E 

03 

0.13b 

04 

0 . 1 7 E 

04 

WAV 

2200.0 

0.0 

0 . 7 5E' 

-03 

0.115 

01 

0.32E 

01 

0.56E 

01 

0.93E 

01 

0.11E 

02 

WAV 

2200.0 

0.25E 

01 

0 . 1 5E 

02 

0.165. 

03 

0.35E 

03 

0.55E 

03 

0.96E 

03 

0. 135 

04 

WAV 

2225.0 

0.0 

0.  1 8 E 

-01 

0 . 32E 

01 

0.66b 

01 

0.93E 

.01 

0 . 1 2 E 

02 

0.12E 

02 

WAV 

2225.0 

0.25E 

01 

0. 1 4E 

0 2 

0.142 

03 

0.25b 

03 

0.40E 

03 

0.7  3E 

03 

0.90E 

03 

WAV 

2250.0 

0. 78E- 

•04 

0.28b 

00 

0.785 

01 

0.12b 

0 2 

0.14E 

02 

0 . 1 4 K 

0 2 

0 . 1 2 E 

02 

WAV 

2250.0 

0,255 

01 

0.14E 

0 2 

0.115 

03 

0 . 1 7E 

03 

0 . 2 ft  E 

03 

0,46b 

03 

0,63  b 

03 

WAV 

2-27  5.0 

0.215- 

•01 

0 . 2 ft  E 

01 

0 . 1 5E 

0 2 

0.1  7t 

02 

0.17b 

02 

0.15E 

02 

0.12  E 

02 

WAV 

2275.0 

0 . 2 5 E 

01 

0. 12E 

02 

0 . 7 5 2 

0 2 

0 . 1 4 E 

03 

0.  20E 

03 

0.37E 

03 

0 . 4 1 £ 

03 

WAV 

2300.0 

0 . 1 9 E 

01 

0.1  6E 

02 

0.2.22 

02 

0.20  £ 

02 

0. 19E 

02 

0 . 1 6 E 

02 

0. 1 3E 

02 

WAV 

2300.0 

0.23  E 

01 

0.78E 

01 

0.852 

02 

0 . 1 1 E 

03 

0.22E 

03 

0. 34E 

03 

0.415 

03 

WAV 

2325.0 

0.37E 

02 

0.36E 

0 2 

0 . 2 3 E 

02 

0.2  OK 

02 

0 . 1 8 E 

02 

0.13E 

02 

0.93E 

01 

WAV 

2325.0 

0 . 1 4 E 

01 

0 . 5 7 E 

01 

0.425 

02 

0.73E 

02 

0.79E 

02 

0 . 1 1 E 

0 3 

0.135 

03 

WAV 

2350.0 

0 . 1 1 E 

0 2 

0 . 2 1 E 

0 2 

0 . 2 3 E 

02 

0.20b 

02 

0 . 1 7 E 

02 

0 . 1 1 E 

0 2 

0.7  2E 

01 

WAV 

2350.0 

0 . 2 8 E 

01 

0.80E 

01 

0.242 

02 

0.34E 

02 

0.37E 

02 

0.49E 

02 

0.56E 

02 

WAV 

2375.0 

0 . 2 4 E 

02 

0.34E 

02 

0.205 

02 

0. 14E 

02 

o . 1 1 e 

02 

0. 60E 

01 

0.27E 

01 

WAV 

2375.0 

0 . 1 1 £ 

01 

0 . 2 1 E 

01 

0.55  2 

01 

0.75b 

01 

0.70E 

01 

0 . 6 8 E 

01 

0.63E 

01 

WAV 

2400.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2400.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2425.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2 A 25 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2450.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2450.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2475.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2476.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2500.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

WAV 

2500.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0  ' 

APPENDIX  III 


SAMPLE  CASE  OUTPUT 


WAVE  NUMBER 

WAVE  LENGTH 

radiation 

WAVE  NUMBER 

WAVE  LENGTH 

RADIATION 

vJ/SR/OI/MICRQN 

2000. 

5.0003 

0.98137 £-03 

2025. 

4.9333 

0. 201 15E-02 

2050. 

4.3730 

0.39446E-02 

2075. 

4.8193 

0.76348E-02 

2 10  0. 

a. 751 9 

U.17107E-01 

2125. 

4.7059 

0. 29450E-O1 

2150, 

4.6512 

0 . 4969o£-0  1 

2175. 

4.5977  . 

0.77929E-01 

2 2 0 0. 

4.5455 

0.11446E  00 

2225. 

4.4944 

0.15660E  00 

2250. 

4.4444 

0 . 1 9604E  00 

2275. 

4.3956 

0.22414E  00 

2 3 0 0 . 

4.3473 

0.25091E  00 

2325. 

4.3011 

0.22054E  00 

2350. 

4.2553 

0.1 /91oE  00 

2375. 

4.2105 

0 . 7 S 1 6 6E-0 1 

2 *1 0 0 . 

4.1567 

0.0 

2425. 

4.1237 

0.0 

2 4 50. 

4.0816 

0.0 

2475. 

4.0404 

0.0 

2500. 

4.0000 

0.0 

total  radiation  emitted 

=0. 8259E 

-01 

WATTS/STER 

«AVE  NUMBER 

WAVE  L E Al 

GTH 

Radiation' 

WAVE  NUMBER 

WAVE  LENGTH 

RADIATION 

2000. 

5.0000 

0.116055.-02 

2025. 

4.9383 

0.23865E-02 

. 2050. 

4.9780 

0.45994E-02 

2075. 

4.8193 

0. 91333E-02 

210  0. 

4.7619 

0.19295001 

2125. 

4.7059 

0.3301 4E-0 1 

2150. 

4.6512 

0.55326E-01 

217  5. 

4.5977 

0.35780E-01 

2200. 

4.5455 

0.12457E  00 

2225. 

4.4944 

0.16815E  00 

2 2 50. 

4.4444 

0.20694E  00 

2275. 

4.3956 

0.23378E  00 

2 3 0 0 . 

4.3478 

0.2S03oE  00 

2325. 

4.3011 

0.22397E  00 

2350. 

4.2553 

0.18512E  00 

2375. 

4.2105 

0.32563E-01 

2400. 

4.1667 

0.0 

2425. 

4.1237 

0.0' 

2-450  . 

4.0316 

0.0 

2475. 

4.0404 

0.0 

2500. 

4.0000 

0.0 

total  radiation  for 

STATION 

4 

= 0.76207E-01WATTS/SR/CM 

PR1-0X 

STATION  NO.=  2 8 AXIAL.  DIST.  CCH)=  .11.035  ASPECT  ANGLE  (DEG)=  90.0 


SOM  UP  SHmDSR  AT  X STATION  0.0 


WAVE  NUMBER 

.4  A ‘2  E L z#  N GT  H 

RADIATION 

WAVE  NUMBER 

WAVE  LENGTH 

RADI  A HON 

W/SH/CM/MICRON 

2000. 

5.0000 

0.0 

2025. 

4.9383 

0.0 

2 0 S 0 . 

4.3780 

0.0 

2075. 

4.8193 

0.0 

2100. 

4.7419 

0.0 

2125. 

4.7059 

0.0 

2150. 

4.6512 

0.0' 

'2175. 

4.5977 

0.0 

2200. 

4.5155 

0.0 

2225. 

4.4944 

0.3 

2250. 

4.4444  ■ 

0.0' 

2275. 

4.3956 

0.0 

2300. 

4.3473  ' 

0.0 

2325. 

4.3011 

0.0 

2350. 

4.2553 

0.0 

■2375. 

4.2105 

0.0 

2400. 

4.1567 

■ 0.0 

2 425. 

4.1237 

0.0 

2450. 

4.0316 

0.0  . 

2475. 

4.0404 

0.0 

2500. 

4.0000 

0.0 

• 

TOTAL,  RADIATION  EF.ITTEO 

=0. 1653E  01 

WATTS/STER 

% 

•tAVE  NUMBEk 

« A V E GEN 3T H 

RADIATION 

WAVE  NUMBER 

WAVE  LENGTH 

RADIATION 

2000. 

5.0000 

0.1 0465E-0 1 

2025. 

4.9383 

0. 21h07E-01 

2050. 

4.3730 

0. 4194 2 E- 01 

2075. 

4.3193 

0. 31130 E- 01 

2100. 

4 . 7o  1 9 

0.26070E  00 

2125. 

4.7059 

0.47098E  00 

2150. 

4.6512 

0.328 9 IE  00 

2175. 

4.5977 

0.13972E  01 

2200. 

4.5455 

0.21785E  01 

2225. 

4.4944 

0.31900E  01 

2250. 

4.4444 

0.42630E  01 

2275. 

4.3956 

0.50712E  01 

2300. 

4.3473 

0.56v53c01 

2325. 

4.3011 

0. 51826 E 01 

2 3 5 0. 

4 . 2 3 5 3 

0 . 3 9 99oS  01 

2375. 

4.2105 

0.16354E  01 

2 4 00. 

4.1657 

0.0 

2425. 

4.1237 

0.0 

2450. 

4.0515 

0.0 

2475. 

4.0404 

0.0 

2 3 0 0. 

4. OOoO 

0.0 

TOTAL,  KAOIA  TIOM  FOR  STATION  28  = 0.0 


wATTS/SR/CM 


I 

I 


3>0  JftM  ~7Q 


LAPP  <$AS  DYNAMICS  — A<W  fcftDtfcTtb^  cod£ 


*K[  - , 


PRJ .ox 


W Isdfr 


o.o 

0.41 
0.61 
1.22 
1 . f>l 

1.57 
2.37 
2.73 
3.18 

3.58 
B . S8 
4.35 
4.32 
5.22 
5.63 
6.03 

1 6.44 

6.85 
6.94 
7.41 

7.85 
8.30 
8.76 
9.21 
9.66 

10.12 

10.58 

11.04 


0.5837254E-01 
0.64682846-01 
0. 7061S87E-01 
0.78207146-01 
0.  d6495286-01 
0. 93884896-01 
0.1016850E  00 
0. 1099497E  00 
0.11H5646E  00 
0.1266 6 3 8E  00 
0.1350876E  00 
0.143d736E  OO 
0.15242866  00 
0. 1602380 B 00 
0.1o65467E  00 
0. 17201  12E  00 
0.17827736  00 
0.18491186  00 
0. 1 851264E  00 
0.19225106  00 
0.19836295  00 
0.20381165  00 
0.2086537E  00 
0.2134077E  00 
0.2167o25E  00 
0.2199754E  00 
0. 2249740 E 0 0 
0.0 


STATION  RADIATION  Cw/SR) 

Pr-4 

0.25220106-01  . 

0.27304626-01 

0 . 3 0 0 6 3 1 2 £ -0 1 

0.32622936-01 

0.32326526-01 

0.391596.56-01 

0.42525826-01 

0.46564486-01 

0.48428206-01 

0.52257146-01 

0.58043456-01 

0.6281853E-01 
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